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Abstract — This paper considers maximizing tlie network-wide 
minimum supported rate in tlie downlinli of a two-ceil system, 
where each base station (BS) is endowed with multiple antennas. 
This is done for different levels of cell cooperation. At one 
extreme, we consider single cell processing where the BS is 
oblivious to the interference it is creating at the other cell. 
At the other extreme, we consider full cooperative macroscopic 
beamforming. In between, we consider coordinated beamforming, 
which takes account of inter-cell interference, but does not 
require full cooperation between the BSs. We combine elements 
of Lagrangian duality and large system analysis to obtain limiting 
SINRs and bit-rates, allowing comparison between the considered 
schemes. The main contributions of the paper are theorems 
which provide concise formulas for optimal transmit power, 
beamforming vectors, and achieved signal to interference and 
noise ratio (SINR) for the considered schemes. The formulas 
obtained are valid for the limit in which the number of users per 
cell, K, and the number of antennas per base station, A^, tend to 
infinity, with fixed ratio /3 — K/N. These theorems also provide 
expressions for the effective bandwidths occupied by users, and 
the effective interference caused in the adjacent cell, which allow 
direct comparisons between the considered schemes. 

Index Terms — linear precoding, regularized zero forcing, inter- 
cell interference, base station cooperation, multicell processing, 
cellular systems, MIMO broadcast channel, interference channel 



I. Introduction 

A. Problem scope 

Consider the downlink (DL) of a cellular system in which 
a base station (BS) services many mobiles within the cell. If 
the BS is equipped with multiple antennas we have the classic 
MIMO broadcast channel (BC), which has been the focus of 
much attention in the past few years, including the recent, 
celebrated characterization of the capacity region using dirty 
paper coding There has also been substantial interest in 
devising suboptimal, but practical approaches based on linear 
precoding techniques (i.e. beamforming). 

The MIMO BC is the appropriate model for a single isolated 
cell. What happens when we bring several cells together, so 
that each is affected by the interference from the others? 
This paper analyzes the performance of linear precoding in 
a multicell setting. 

Consideration of interference leads us to examine the system 
level architectural issue of cooperation between nodes in the 
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network. The traditional approach to interference is to partition 
the cells in time or bandwidth to avoid a strong interference 
coupling between adjacent cells. However, with multiple an- 
tennas at each base station, this may be suboptimal. When 
multiple antennas are incorporated at each BS, we can tradeoff 
the maximization of the rates of the in-cell users (ignoring 
interference), with the minimization of the interference spilled 
over into the other cells. If enough cooperation is enabled, 
these two objectives can go hand in hand. 

A large body of research has recently dealt with coopera- 
tion and coordination in multicell systems. Many papers are 
concerned with developing new algorithms, to meet various 
proposed performance objectives (e.g. transmit power mini- 
mization under given SINR constraints for coordinated beam- 
forming (CBf) Q, I3J, or minimum SINR maximization for 
network MIMO fA\). Others provide a performance analysis 
of a given scheme under a particular channel model. 

One can distinguish between scenarios where BSs each 
serve a different group of users, and cases where all the 
transmitters jointly transmit to all users in the system, the so- 
called multicell processing (MCP), macrodiversity or network 
MIMO. For MCP, a classical model for performance analysis 
is Wyner's model ISj, ||6]. This model was first used on the 
DL in Q, where a linear pre-processing dirty paper coding 
approach is proposed. In ||6l, the sum rate is characterized 
for the case where single-antenna base stations pool their 
antennas together to perform zero-forcing (ZF) to the users in 
the system. More precisely, they consider a circular variant of 
the infinite linear Wyner model for both non-fading and fading 
scenarios, with scheduling based on local channel statistics. 
Results are derived for the regime in which the number of 
cooperating BSs tends to infinity, and scaling results are also 
obtained by letting the number of users per cell do the same. 

In this work, we focus on optimizing linear precoding under 
different states of CSI and data sharing. We assume that each 
BS can simultaneously serve more than one user, and we 
formulate the problem of maximizing the minimum network- 
wide achievable rate, i.e. rate balancing, under the following 
three different architectures: 

i) Single cell processing (SCP), in which each BS has 
perfect CSI about mobiles in its cell, but no knowledge 
about the channels to mobiles in other cells. Here, the 
BS can control interference between the mobiles in its 
own cell, but not the interference spill-over to other cells. 

ii) CBf, in which the beamforming decisions at each BS 
take into account the impact of interference on the other 
cells. In this approach, each BS has CSI about the 
channels to its own users, but it also knows the channels 
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to the mobiles in the other cells. This allows the BS to 
control the interference it causes to the mobiles in other 
cells. 

iii) MCP, in which the BSs cooperate to jointly precode 
signals to all the mobiles in all cells. Through coop- 
eration, the BSs know the channels from all BSs to all 
the mobiles in the whole network, and hence they can 
jointly beamform, using their combined resources. 

In the first two approaches, user data is routed to a single BS 
only, whereas in the MCP case, it is routed to all cooperating 
BSs. We specialize our derivations to only two cells. Never- 
theless, such an approach can be extended to larger networks, 
if one wishes to perform system-wide optimization of a large 
cellular network. This is the object of ongoing research f8|. 

Even the two cell model gets very complicated when we 
have multiple antennas at the BSs and many mobiles in each 
cell. With independent flat fading between each transmitter- 
receiver pair, many parameters need to be specified. For this 
reason, we assume that they are selected randomly, and we 
take a large system approach in which the number of antennas 
at the base station, N, and the number of mobiles in each cell, 
K, both grow large together, while the ratio ^ , which we refer 
to as cell loading, is held at a constant, denoted /3. 

B. Related work 

Random matrix theory has received a lot of attention in 
the communications literature recently |[9l, particularly large 
system results where some of the system parameters, such 
as number of users, the length of the signature sequence (in 
CDMA scenarios), or the number of transmit and receive 
antennas (in MIMO settings), are allowed to grow large at 
the same rate. Such asymptotics often produce a compact 
characterization of performance in the large system limit, 
amenable to system optimization. An attractive feature of 
these results is that the asymptotic expressions turn out to be 
good estimates for even relatively low values of the scaling 
parameters flo!, HI, lUl, lUl, M . 

Most applications of large system analysis have been for 
the uplink (UL). Surprisingly, there has been limited work 
on the DL until quite recently, despite a well established UL- 
DL duality theory. Recently, however, a large systems analysis 
of regularized ZF (RZF) beamforming was carried out to 
characterize its limiting performance in a single cell context, 
allowing the optimization of the regularization parameter ifTSl . 
and II6 1 considers ZF and RZF for correlated channel models. 
The present paper generalizes 1 15 1 to the multicell context, and 
to a wider class of beamformers, exploiting an UL-DL duality 
theory. 

In the past year or so, there has been further work that 
explores the interplay between UL-DL duality, Lagrangian 
optimization and large systems analysis. In |17|, duality be- 
tween MAC and BC is used to characterize and optimize 
asymptotic ergodic capacity for correlated channels. A similar 
approach is taken in IITSl . llT9l . EOll to treat a large systems 
analysis of ergodic, weighted sum-rate maximization. These 
papers consider maximizing network utility to obtain fairness, 
under user scheduling, and they consider various forms of 



base station cooperation (clustering). The aim is to use large 
system analysis to obtain numerical methods that are much 
more efficient than Monte-Carlo simulations; in ll20l they 
obtain an "almost closed-form" numerical analysis tool. Linear 
ZF beamforming and non-ideal CSIT are considered in flQl. 
In 121] random matrix theory is applied to a different CBf 
setup than the one considered here: more particularly, they 
consider the problem of weighted sum of the transmit powers 
minimization CBf problem initially formulated in |2|, and 
propose a strategy which also requires instantaneous local CSI 
and sharing channel statistics. Monte Carlo simulations are 
resorted to in order to claim asymptotical optimality of the 
results; these are however derived for a more general channel 
model than we use in the present paper In f22|, we provided 
preliminary results that we now present in more detail. 

There has also been some very recent interest in using 
large systems analysis to succinctly answer questions con- 
cerning channel uncertainty, optimal amount of training, and 
required rates of feedback of channel state information. These 
questions are addressed for the downlink of a single cell in 
II23I . where deterministic equivalents of SINRs for ZF and 
RZF are derived, under channel uncertainty, and the resulting 
expressions are then optimized with respect to the number of 
users and the number of symbols devoted to training. In |24|, 
large systems analysis is used to optimize the number of BSs 
that should be cooperating on the UL, taking into account 
unreliable backhaul links, and the channel estimation required 
to measure the channels. Both papers consider much more 
general channel models than we do in this work. 

Many duality results have been established in the context of 
MIMO communications. The first UL-DL duality result was 
obtained for the point to point MIMO channel in 125). Another 
early work, |26|, considered joint optimization of power and 
beamforming vectors for the DL of a multiple antenna cel- 
lular system employing a simple linear transmission strategy 
followed by single user receivers, such that the SINR at each 
mobile is above a target value; the proposed algorithm achieves 
a feasible solution for the DL if there is one and minimizes the 
total transmit power in the network. In ||26| . connections were 
made between this problem, and the UL power control frame- 
work of 127], and we exploit these connections in the present 
paper (as have many other authors). In |28|, the problem of 
DL power minimization subject to target SINR constraints 
is addressed using a Lagrangian duality framework, and the 
transmit powers of the dual UL are found to correspond to the 
Lagrange coefficients associated with SINR constraints. 

In the context of the capacity region of the Gaussian BC, 
duality results are provided in |291, JSOl, EB, EH. When 
linear beamforming is considered, as in the present paper, 
several duality results have been obtained, and applied to 
designing iterative algorithms for DL beamforming. It was 
shown in |33| that, under a sum power constraint, UL and 
DL achievable rate regions are the same. Effective bandwidth 
results derived for the UL were also proven to hold in the DL. 
A similar approach to duality for the BC is taken in If34l . 

Over the past decade, a large body of research has appeared 
in which beamformers for the BC are derived via Lagrangian 
optimization techniques. Optimization formulations include 
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the minimization of the sum of mean-squared errors ll35l . 
|l36l, fy7\, f3S\, the minimization of power subject to SINR 
constraints |39|, and the maximization of SINR's subject 
to power constraints 1401 . |f39l . In 11411 . the problems of 
maximizing the sum of effective bandwidths in both the UL 
and DL (via duality) are considered. For the two-user case, 
a closed form solution may be obtained, whereas the more 
general problem can be formulated as a convex optimization 
problem and solved via an interior point method. Related 
problems such as SINR balancing are also considered. In PBl . 
the Lagrangian duality approach is extended to include per- 
antenna power constraints. It is shown that the dual problem 
corresponds to an UL with uncertain noise at the receiver, and 
the focus is the derivation of various algorithms for solving 
the power minimization problem. 

Some of the above-mentioned formulations are not directly 
convex problems; for example the DL beamforming problem 
to minimize power subject to SINR constraints is not a convex 
problem. However, it can be transformed into a convex prob- 
lem, as shown in f39l, which shows that strong duality holds 
in the original non-convex problem [42 J. Thus Lagrangian 
methods can be applied to these non-convex problems, and 
this fact is exploited in the present work. 

The optimization approach to beamforming provides much 
insight into questions of system capacity, but when it comes 
to practical beamforming design other considerations also 
come into play. For example, using an iterative algorithm to 
find an optimal beamformer may not be viable if channels 
change rapidly over time|^For this reason, many authors have 
focused on beamformers that can be found without iterative 
methods, such as the classical ZF beamformer]^ Of course, 
rather than adapt to the instantaneous channel realizations, an 
alternative approach is to adapt to changing channel statistics. 
Stochastic approximation methods have a long history, and 
provide iterative algorithms for these scenarios. It is of interest 
to extend the results of the present paper to smaller (non- 
asymptotic) systems and such methods may prove useful. 

ZF beamforming steers nulls at the other users so that they 
each get an interference-free version of their desired signal. 
However, the signal to noise penalty can be very high if the 
matrix to invert is ill-conditioned. RZF beamforming adds a 
regularization term to the ZF beamformer to provide numerical 
stability and better performance |44|. As noted earlier, |[T5l 
undertakes a large system analysis of RZF, exploiting how its 
beamformer resembles the LMMSE receiver on the UL. 

A number of works have recognized that the optimal beam- 
former can have a RZF structure in some special cases |35|, 
(361, 1371, EOl, im, 1391 ■ Note that for the power minimiza- 
tion subject to SINR constraints problem |39|, the simple form 
of the regularized beamformer only arises when there is a great 
deal of symmetry in the channel model. The examples given 
in Section VI of ll39ll . where such a form emerges, consist 

'The classical Foschini-Miljanic|'43 1 power control algorithm attempts to 
find the optimal power levels in a dynamic environment, but with beamforming 
there are many more parameters to be determined. 

^By "iterative" we mean schemes that require iterative updates of physical 
parameters, such as transmit powers and beamforming vectors, not simply 
numerical methods that can be implemented within a single network node, 
e.g. in the computation of a matrix inverse. 



of the diagonal case (i.e. no interference between channels) 
and the symmetric case, when the channel matrix has equal 
diagonal elements, and equal off-diagonal elements. In the 
present paper, channel matrices are randomly selected, so this 
particular symmetry condition is not satisfied. However, our 
model is symmetric in the average statistics of the channel 
matrices, which leads to the RZF beamformer resulting from 
the large systems limit. Nevertheless, this form of symmetry 
is more general than the very special one considered in |39|. 

In general, ll39l uses iterative techniques to find numerical 
solutions for the optimal beamformer using conic optimization 
techniques. The latter are exploited in |42| to handle per- 
antenna power constraints, and extended further in |[2l, Q 
to provide a duality theory for multicell systems in which 
there are per-BS power constraints. In Q, lO, CBf, a novel 
beamforming strategy in which the BS takes into account the 
interference it creates in adjacent cells, is proposed. f3l focuses 
on developing fast algorithms to find optimal beamformers, 
and numerical evidence is provided to show the improved 
performance due to coordination, relative to traditional SCR 

The RZF structure emerging from our symmetric two cell 
model is reminiscent of beamformers designed to optimize 
other criteria, such as minimum variance distortionless re- 
sponse (MVDR) beamformers. In particular, there is a sim- 
ilarity between the RZF structure and the so-called diagonal 
loading structure that has been used to deal with uncertainty 
in the noise correlation matrix |45|. In |46|, a large systems 
analysis of MVDR beamformers with diagonal loading is 
undertaken where the number of antennas, and the number of 
observations, grow large together. Deterministic equivalents 
are obtained, allowing the diagonal loading factor to be 
optimized. 

In Section [III] we review the optimization problems asso- 
ciated with the classical (non-coordinated) SCP, and the CBf 
strategy of |3|. We also formulate the optimization problem 
for MCP in which the BS antennas jointly coordinate their 
transmissions, as in network MIMO |47|. 

The main contribution of the paper is a large systems 
analysis of these three schemes, and optimizations based on 
the limiting asymptotics. In particular, we provide a clean 
characterization of each of the three schemes in the large 
systems limit in terms of effective bandwidths. In the case 
of SCP, we show that the RZF beamformer is asymptotically 
optimal in the large systems limitj^ In the case of CBf, we 
show that the asymptotic form of the optimal beamformer 
provides an enhanced version of regularized ZF, one that 
we call generalized regularized zero-forcing (GRZF)0 This 
beamformer is a novel contribution of the present paper 
Notation: Lowercase letters, boldface lowercase letters and 
boldface uppercase letters are used to represent scalars, vectors 
and matrices, respectively, x^and are the transpose and 
conjugate transpose of vector x, respectively. CAf{fJ,, C) de- 
notes a circularly symmetric complex Gaussian random vector 
of mean fj, and covariance matrix C. Finally S ^ means S 
is a positive semidefinite (PSD) matrix. 

^^the structure that emerges is more general than RZF, since it also applies 
in regimes in which ZF is not defined. 

■* again, it also applies when ZF is not defined. 
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II. System Model 

Our model, illustrated in Figure [T] has two cells, and each 
BS is equipped with an array of N antennas. There are K 
single-antenna mobiles in each cell. We assume flat fading and 
adopt the following notation regarding channel coefficients: 

i) the channel vector from BS j to user k in cell j' is 



denoted hk.j'j, where hk^ 



ii) the channel vector from all BSs to user k in cell j' is 
denoted hk.j'', in other words, hkj' = [hfc./,i hj. j' 2]- 

The data symbols to be transmitted to each user are assumed 
to be independent identically distributed (i.i.d.) CAf{0, 1) ran- 
dom variables (r.v.), and the data symbol for user k in cell j 



is denoted 



Let s„ 



5ij 



SKj]'^ and s — [sf sj]'^ 



The received signal at user k in cell j is given by 

2 



J2 ^^^jj" 



(1) 



where Xji G C denotes BS /'s transmit signal, which 
consists of the linearly precoded symbols of the users it is 
serving, and which is subject to the average power constraint 
E [x|fxj'] < P and nk,j ^ CJ\f{0,a^) is the receiver noise. 
Receivers perform single-user detection, i.e. treat interfering 
signals as noise. The way the precoding vector Xji is generated 



depends on the considered scheme; See Section III 



We assume the channels between each BS and user are in- 
dependent. Moreover, channels between a user and his serving 
BS are i.i.d. CAf (0,1) whereas channels between a user and 
an interfering BS are i.i.d. CJ\f{0, el). Thus, e controls the 
interference level between neighbor cells. This is a simplified 
model of a cellular network, yet provides useful insights 
It is the two cell DL version of Wyner's model P9ll . 




Fig. 1. System model 



III. Cooperation and Linear Beamforming 
Schemes: Primal problems 

We consider the problem of maximizing the network-wide 
minimum achievable rate for three different degrees of cooper- 
ation and coordination between the cells. The optimization to 
be carried out in each case is presented. The first problem is the 
now classic one treated in |26|, the second is the cooperative 
scheme proposed in [i3j, and the third is new, although similar 
approaches have been proposed in several papers (e.g. |47|). 
In the equations below, when we consider a particular BS j, 
index j = mod {j, 2) + 1 corresponds to the other BS. 



A. Single cell processing 

This is the conventional case where each BS serves its own 
users without worrying about the other cell. Here, we assume 
full re-use of time and spectrum across cells. Xj is of the form: 



K 

E 

fe=i 



(2) 



where the symbol for user k in cell j, Skj, is linearly precoded 
by beamforming vector w^j. is the concatenation of the 
K precoding vectors in cell j into a N x K matrix, the fcth 
column being Wkj- 

In cell j the problem to be solved is the following: 



max. 



7,wfcj,fe=i,...,_fs: 7 



s.t. SINRfc,j > 7 

K 



fe=i 



wfejir < p. 



The SINR at user k in cell j is given by 

|hfe,ijWfej|2 



SINR 



(3) 



(4) 



al is the noise plus other-cell interference power at that user: 



K 

k' = \ 



(5) 



and needs to be fed back to BS j. 

Solving Q may require an iterative procedure, since beam- 
forming at each BS influences the interference, and therefore 
the transmission design, at the other. We focus on finding the 
maximum SINR that can be met in both cells. This can be 
obtained using a bisection method. Thus, for fixed 7 we obtain 
the beamforming vectors by minimizing total transmit power 
subject to SINR constraints on the cell's users. If the optimum 
is < P for both BSs, 7 is attainable. We thus focus on solving: 

K 



mm.w,,j=i,...,K ^\\^kj\\- 

k=\ 

S.t. SINRfej > 7, k^\,...,K. 



(6) 



For 7 achievable with unlimited transmit power, its solution 
will be a set of beamforming vectors that minimizes the 
total cell transmit power, and achieves SINR 7. Further, in 
Theorem [T] we will provide a condition that is both necessary 
and sufficient for the target SINR to be achievable, given 
unlimited power The maximum 7 is the target SINR for which 
the transmit power constraint is met with equality. 



B. Coordinated Beamforming 

Here, each BS sends data to its own users only, as in SCP, 
but CSI is shared between the BSs so that the interference 
generated in other cells is taken into consideration. Xj will be 
similar to that in (|2]i, although the precoding design differs. 
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The problem formulation in (|3]l becomes 

max.^^wfcj,fc=i,...,i<-,i=i,2 7 

s.t. SINRfe,,>7, fc = l,...,i^,j = l,2 

K 

Y.\\^kjf<P, J = 1,2, (7) 

k=l 

which is a joint, two-cell optimization problem, requiring a 
coordinated solution. The SINR at user k in cell j is similar 
to Q but can be expanded into 

|hfcjjWfcj|2 



SINR 



k.j 



(8) 



since all the channels are centrally known. 

Here too the problem may be solved by a bisection method. 
To determine feasibility of a given 7, we solve (as in LSjJ^ 

min.0>o,wfcj,fc=i,...,_R',j=i,2 2P0 

s.t. SINRfcj>7, fc = l,...,X,j = l,2 



is' 



(9) 



where (f)P upper bounds the power expended at each BS, 
and the aim is to use the minimum power level. Clearly, the 
maximum network-wide achievable SINR is the 7 for which 
the optimal (j) — 1. 

C. Multicell processing 

This is the case where the BSs cooperate fully: both CSI and 
data is available at all transmitters, who pool their antennas 
together to serve the users jointly. The transmitted signal x — 
[xi; X2] will be of the form 



2 K 

x-^^WfejSfcj = Ws, 

3=1 fc=i 



(10) 



where W e (^2^x2^ j-j^g overall precoding matrix. The 
MCP optimization problem is 



s.t. SINRfej>7, fc = l,...,i^,j = l,2 

2 K 

^ ||E,Wfc,,,.|l'<^, J = l,2, (11) 
i'=ife'=i 

where matrix Ej , j = 1,2 is diagonal and used to select the 
elements of each beamforming vector corresponding to BS 
j (i.e. its non-zero diagonal elements occupy locations (j — 
l)iV + 1 to jN). The MCP SINR at user k in cell j is 

SINR,, = . . (12) 

^Note that the constant 2P factor in the objective function is included as 
it leads to P being eliminated from the dual problem formulation. 



Once again, the above problem may be solved by the 
bisection method. To determine feasibility of a given 7, we 
solve, as in the CBf case, the following optimization problem; 

min. 2P4> 
s.t. SINRfcj>7, fc l,...,i^,j = 1,2 

2 K 

J2 E w^j^k'j'f <^p, j = i.2. (13) 
j'=ifc'=i 



IV. Solution VIA Duality THEORY 

The above optimization problems are non-convex, and ap- 
parently nontrivial to solve. Similar problems of power control 
and receiver optimization on the UL are easier to solve, since 
each UL receive vector can be individually optimized. As 
discussed in Section |I-B| motivated by this observation, a 
number of duality results have been established connecting 
DL optimization problems, to corresponding dual UL ones. 
Most relevant to us is the elegant duality theory developed 
by Yu and Lan ll42l to handle per-antenna power constraints, 
exploiting earlier work on conic optimization 1391 . 

Algorithms analogous to those proposed in p2 l can be 
applied to solve the above problems, but even after these 
iterative algorithms have converged to the optimal solutions 
there is a great deal of complexity in the form of the optimal 
beamformers. Our goal is to use the duality theory to derive 
suboptimal (but, in the context of our simplified model, asymp- 
totically optimal) beamformers that have relatively simple 
structures. In fact, generalized regularized beamformers will 
emerge from our analysis. Our approach will be to apply a 
large systems analysis to the dual, virtual UL optimization 
problems. Following ll42l . we begin by writing down the dual 
problems to (|6]l, (|9]l, and ( pjj ), respectively. 



A. Dual UL SCP problem 

Let Xkj/N > denote the Lagrange multipliers correspond- 
ing to the SINR constraints in (|6|. The dual problem is 



K 



k=l 



k'^k 



k = l,...,K. (14) 



As explained in |42|, strong duality holds despite non- 
convexit}]^ so the optimal value in the dual problem is equal 
to that of the primal DL beamforming problem. 

The dual variables Xkj/N may be thought of as dual UL 
transmit powers, in the following dual UL power control and 

^because the non-convex problem can be transformed into a convex problem 
using the techniques in 1391 
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beamforming problem ll42l : 



mm.w^^,Afcj,fc=i,...,if ^ 



s.t. max ■ 

Wfcj 



>7, 



fc = l,...,if. (15) 

Indeed, ([14]) can be shown to have the same optimal values as 
^T5[ . The optimal w^j, up to a scalar, are given by 



E 



(16) 



which we recognize to be MMSE filters for the UL problem. 

To see why these problems are equivalent, note that the 
SINR constraints in the UL problem ( fTS] ) become, after 
substituting for the MMSE filters: 



and the {pkj)^^i can be determined from this set of equations. 
We undertake a large system analysis of this scheme in 
Section |V] and present the results in Theorem [T] 

As a final remark, we note that the UL-DL duality described 
in this section is not the same as the celebrated UL-DL duality 
used to characterize the capacity region of the MIMO BC. In 
the present section (and indeed throughout this paper), we do 
not allow so-called dirty paper coding (DPC), which is the 
coding technique that achieves the capacity of the MIMO BC. 
The original notion of UL-DL duality was a correspondence 
between points in the capacity region of the MIMO BC with 
points in a dual UL multiple access channel with suitably 
chosen noise covariances. However, it has been shown that the 
optimization approach developed in f42\ to obtain the optimal 
beamforming vectors for the DL power minimization problem 
(as reviewed in the present section) can be extended to allow 
DPC, and hence to obtain the MIMO BC capacity results, 
including the setup with per antenna power constraints (see 
iS). 



TV 



J2 -j[-^k\jj^k',j.j 



k'^k 



h^«>7. (17) 



The optimal AfejS are the unique solution to the fixed point 
equation: 



(18) 



Since 7 is achieved exactly, the SINR constraints in ([TSj are 
achieved with equality. It follows that if the minimization in 
([TSj is changed to a maximization, and the SINR inequalities 
are reversed, then the same solution will be found. It is shown 
in B2I that in so doing, one obtains ( [T4| l. 

The right hand side of ([T8| is an interference function of 
the dual UL powers, in the sense of Yates' framework on 
UL power control ll27l . Thus, iterative approaches to power 
control converge to the optimal dual variables. Indeed, this was 
the original approach to the multicell beamforming problem 
taken in ll26l . one of the first papers to exploit an UL-DL 
duality. More recently, more efficient approaches to solving the 
dual problem have been proposed [42|, but these still require 
several iterations before the algorithm converges, which could 
be problematic. The large systems analysis in Section[V] allows 
us to derive simple yet asymptotically optimal beamformers. 

From the Karush-Kuhn-Tucker (KKT) conditions, one can 
show that the optimal w^j for the primal are the same, up to 
a scaling factor, as w^j. Thus, they can be written as 




-j^ is the transmit power allocated to beamforming 
vector -Wkj on the DL. From the DL SINR constraints, for 
k = l,...,K, 



Pkj |Jl/c,j,jWfcj| 



k'^k 



N 



\vfk'j\ 



'k,j^ 



(20) 



B. Dual UL CBf problem 

Letting Xkj/N > denote the Lagrange multipliers corre- 
sponding to the SINR constraints and /ij those corresponding 
to the power constraints in (|9]), its Lagrangian dual problem is 

EAfcj 2 
N 

hk 



max.A,,>o,u,>o 



s.t. /ijl 



\k 



kj ,H 



X 



{k'.rmk,]) 



fc = l,...,if,j = l,2 



2 

E 



(1 - ^lj)P = 0. 



(21) 



As for the SCP case, this problem can be shown to be 
equivalent to a dual UL problem with uncertain noise (cf. 
|42J), or equivalently 

max.^^>omin.A,^.>o,w,, ^'^^ 

k,j 



S.t. 



'kj 



> 



2 

E 



(1 - = 0. 



fc = l,...,X,j = l,2 



(22) 



(19) The Lagrange multipliers /ii,/i2 can be interpreted as noise 



levels at the two BSs, respectively, on the virtual UL B2II . 

For any choice of dual variables, {/ii, ^2), (;^-^fcj)fc-i ^-i' 
the SINR achieved on the virtual UL can be shown to be: 



N 



^k.j,] 



fljl 



E 



A 



k'j' uH 



{k',3')^ik,j) 



N 



(23) 



7 



Using the optimal dual variables, all ( |23| ) e^ual 7, which 
allows us to write down the optimal {^kj)k^ij=i the 
unique solution to a fixed point equation: 



7A^ 



(24) 



H 



(32) 



Once (/ii,/i2) are determined, this equation provides an im- 
plicit solution to the dual problem, in that standard iterative 
methods can be applied to obtain the dual UL powers. It 
remains to find fii and ^2- This is addressed in Section [V] 

The optimal beamforming vectors on the dual UL, w^j, 
assuming feasibility, are (up to a scaling factor) 



E 



k'j' uH 



N 



(25) 



the optim al b eamforming vectors on the DL are of the form 
Wfcj = ^ ^ II II , where is the power allocated to the 
beamforming vector Wkj ■ From DL SINR constraints, for j — 
l,2,k = l,...,K, 



Pkj |hfcj,j-Wfcj| 
^7 l|wfc,l|2 



^ N 

j' ,k' ,{k' ,j')^{kj) 



Pk'j' \hk j,j'-Wk'j'\ ^ ^2 



Wfe' 



(26) 



where hkj = hfc.jM^^/-^, and let hkj is the diagonal matrix 
whose diagonal entries vector is given by 



[A 



; Aa'jS Ai^ ; 



; AkjJ 



(33) 



Using this notation, for any choice of dual variables, {^1,^2), 

K, 2 

, the SINR achieved on the virtual UL (i.e. the 

. k=i,j=i 

left-hand side of Eq. pi] )) can be rewritten as 



N 



N 



(34) 



With the optimal dual variables, all ([34]) ecjual 7, which 
allows us to write down the optimal {^kj)i;^ij^i as the 
solution of a fixed point equation: 

_^ -1 
\H^,Lfe,,Hfe,,+lj hi,] . (35) 



Afe, = h 



kj 



N 



C. Dual UL MCP problem 

With ^'s and /ij's defined as in ( |IV-B| i, the Lagrangian 
dual to problem ( [T3] ) is equivalent to 



max.^,>omm.A,^.>o,wfc, 2^2^ 



3=1 fc=i 



s.t. 



(k',]''mk,3) N "fe'j 



> 



E(i-^j)^o- 



(27) 

7 

(28) 
(29) 



where, to simplify the equations, M = X]j'=i Mj'Ej' is used. 
The optimal w^j is 



^ ^ Afc,,.h|(,^.,hfe,,,. +M 
Plugging this into the inequality in ([28]), we obtain 



(30) 



TV 



E 



(fc'j')#(fcj) 



^ ^■ 



(31) 



Define j as 



As for CBf, once Hi,fJ,2 are obtained, the optimal 
('^*:j)fcLi^j=i implicitly given as the unique solution of 
( [35] ). As the duality gap is zero for a feasible primal, at the 
optimum. 



2P0 : 



^ 2 K 



N ■ 



l^l^^k3^- (36) 

J = l k=l 

Moreover, the beamforming vectors in the original problem 
and in the dual problem are related as follows: 



lPk3 ^ 
|Wfc, 



(37) 



Plugging these into the DL SINR constraints provides the 
solution for the pfc^ : 



Pk3 IhfcjWfcjl 



fk3\ 



E 



Pk'3' |hfcjWfc/jv 



{k',3'WkJ) 



N 



Wfe/ 



(38) 



V. Large system results 



We proceed to a large system analysis of the solutions 
to the above optimization problems, in the limit as TV, K 
grow large, keeping the cell loading ratio P — ^ fixed. The 
basic idea is to apply large systems analysis techniques to 
the dual UL problems, guided by the analysis of similar UL 
problems in the literature ifTOl . Il50l . For example, the SINR 
expression in (VJ) is typically analyzed by considering the 
limit of the empirical distribution of the eigenvalues of the 

matrix (l + Efe'^fc Afe'jhi(^ jh^, j- j) as TV, X ^ 00, 
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with K/N held at /3. Once this is characterized, a law of 
large numbers for the trace of this matrix can be obtained, 
and from that, the limit of the SINR expression in ( [T7| ). 

In the present paper, extra technical difficulties arise. To 
directly apply the known theorems about large random ma- 
trices, we need to assume that the dual UL powers \kj are 
independent of the channel matrix parameters, and are chosen 
independently of each other from some fixed distribution. 
However, being solutions of the UL dual problem, they form 
a collection of dependent random variables which, moreover, 
depend on the channel matrix parameters. Our approach to this 
problem is to provide lower and upper bounds to the SINR's 
using deterministic dual UL powers. We obtain these bounds 
by exploiting the underlying monotonicity structure of the UL 
power control problem |51|, |27|, \26\. With deterministically 
chosen UL powers, we can perform a large system analysis, 
and then provide a sandwich type argument to show that the 
optimal dual UL powers must converge to deterministic values. 
Since this part of the paper is somewhat technical, we relegate 
it to the appendices, and focus on the results of the analysis. 

A. Asymptotically Optimal Beamfonners 

Theorem 1 (Asymptotically optimal beamforming for SCP). 

Assume j3 ^ + 67^ < L Then, asymptotically, SINR 7 is 
achievable at each mobile terminal in the limit as N ^ 00 
with ^ ^ (3 > 0. In this case, the empirical distribution of the 
(normalized) UL dual power levels (i.e. the Xkjs) converges 
weakly to the constant X given by 



A 



7 



1 



(39) 



+7 



the empirical distribution of the ( normalized) DL power per 
users (i.e. the pkjs) converges weakly to the constant 



P 



9 



1-/3 



7 

1+7 



£7 



(40) 



The per BS power converges to P ~ f3p. The asymptotically 
optimal form of the DL beamformer for user k in cell j is 

t,SCP 



,SCP 



'kj 



where wfj^^ 



k'^k 



(41) 



Finally, the asymptotic SINR, 7, is related to the other vari- 
ables via the fixed point equations 

1 1 



7 



1+7 



p ^ 1+7 



(42) 



Conversely, if /3 (^jip^ + £7^ > 1, then, asymptotically, the 
SINR target 7 is not achievable under the SCP strategy. 

Proof See Appendix [P] ■ 

Corollary 1 (SCP). Subject to the per BS power constraint P, 
the maximum asymptotic network-wide achievable SINR for a 



given cell loading factor j3 is the unique positive solution to 
the following fixed point equation: 

1 1 



Iscp = 



(43) 



l+7Scp 

which has the explicit solution "fscp ^?Mfl/ to 



1 



+ e 



1 



Theorem [T] is interesting in that it relates the solution to 
the optimization problem (|6| to a notion of regularized zero- 
forcing proposed in [44 J as a practical approach that is as 
simple to implement as ZF, yet with better performance. It 
was studied in the asymptotic regime in 1 15 |. The beamformers 
defined in ( |4T] i asymptotically lead to a precoding matrix 



W 



SCP 



H 



H 



(44) 



where Cj ensures the power constraint at BSj is met with 
equaUty, and Hjj is the concatenation of channels between 
cell j users and their serving BS. We should add that the opti- 
mal beamformer always exists, even when the ZF beamformer 
does not, so technically we should only refer to it as RZF in 
those scenarios where the ZF beamformer exists. 

Another interesting observation is that Theorem [T] provides 
a condition that is both necessary and sufficient for the target 
SINR, 7, to be achieveable. We can interpret 



7 



as the 

1+7 

effective bandwidth of a user in cell j, and 67 as the effective 
bandwidth of an interferer in cell j. Effective bandwidths 
provide a simple metric by which different beamforming 
schemes can be compared, as shown in the next two theorems. 



Theorem 2 (Asymptotically optimal beamforming for CBf). 
Assume /3 ^ + j^f^ ) < 1- Then, asymptotically, SINR 7 



is achievable at each mobile terminal, in the limit as N ^ 00 
with — >■ /3 > 0. /n this case, the empirical distribution of the 
(normalized) UL dual power levels (i.e. the XkjS) converges 
weakly to the constant X, given by 

7 



A 



^7 
l+e7 



(45) 



the empirical distribution of the ( normalized) DL power per 
users (i.e. the pkjS) converges weakly to the constant 



P 



Xa' 



2 

f7^7 



l+e7 



(46) 



The per BS power converges to P ~ f3p. The asymptotically 
optimal form of the DL beamformer for user k in cell j is 

^Coord 

(47) 



w 



Coord 
kj 




where w^/'"''' 



A 

TV 



E 



^k',j,j^k',3,j 
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Finally, the asymptotic SINR, 7, is related to the other vari- 
ables via 



1 = 



Conversely, if (3 



1 


CP ^ 


P 


1+67 


1+7 




L «7 


V 1+7 ^ 


^ l+e7 



1 



1+C7 



/3 
1+7 



(48) 



> 1, f/jen asymptotically 
the SINR target 7 is not achievable under the coordinated 
beamforming strategy. 

Proof: See Appendix |E] ■ 

Corollary 2 (CBf). Subject to per BS power constraint P, 
the maximum asymptotic network-wide achievable SINR for a 
given cell loading factor /3 is the unique positive solution to 
the following fixed point equation: 



ICBf 



1 



1 



+ 



(49) 



l+7£ 



CBf 



In other words, "fcBf '"^^^ '^f ^ cubic equation. 

Theorem |2] is interesting in that it provides a novel form 
of RZF beamformer|^ This beamformer is not a direct reg- 
ularization of the standard ZF beamformer in a single cell. 
Rather, it is a regularization of a beamformer that zero forces 
the interference it creates at users in the other cell as well 
as its own; in other words, it transmits to a subset of the 
users that it is zero forcing. We call this a generalized RZF 
beamformer, and it is a novel contribution of the present paper 
It can be used in a finite system, where it is suboptimal, 
but relatively straightfoward to implement. Theorem |2] shows 
that it is asymptotically optimal in the class of coordinated 
beamforming strategies. Note also the clean characterization 
of effective bandwidth for this beamformer, and that it provides 
a significant reduction in the effective bandwidth of the other- 
cell users, compared to SCP, when e is non-negligible. 

Theorem 3 (Asymptotically optimal beamforming for MCP). 
Assume P^j^ < 1- Then, asymptotically, SINR 7 is achievable 
at each mobile terminal, in the limit as N ^ 00 with 
^ ^ (3 > 0. In this case, the empirical distribution of the 
(normalized) UL dual power levels (i.e. the Xkjs) converges 
weakly to the constant A given by 



X = 



1 



7 



1 



(50) 



the empirical distribution of the ( normalized) DL power per 
user (i.e. the Pkjs) converges weakly to the constant 



1 



9 

CT 7 



1 



(1-/3^ 



+7 



(51) 



The per BS power converges to P ^ [ip. The asymptotically 
optimal form of the DL beamformer for user k in cell j is 



MCP 



MCP 




MCP\ \ ' 



(52) 



where w^^*^^ 



A 

N 



E 



(53) 



Finally, the asymptotic SINR, 7, is related to the other vari- 
ables via 



7 



1 

(l+e)A ^ 1+7 



1 



(l+e) 



1+7 



(54) 



Conversely, if (3 > 1, then asymptotically the SINR target 
7 is not achievable under the MCP beamforming strategy. 

Proof See Appendix |F] ■ 

Corollary 3 (MCP). Subject to per BS power constraint P, 
the maximum asymptotic network-wide achievable SINR for a 
given cell loading factor (3 is the unique positive solution to 
the following fixed point equation: 

1 1 

7a/cp 



^ (l+e)P ^ l+7l,cp 

In other words, ^*mcp is equal to 



(55) 



(1 + £)P P 



4+1 



(H 



i + 1 



P 



(1 + £)P 

Although MCP is a complex strategy, in that the cooperation 
between BSs is much greater, the reward is better performance 
than that attainable in a single isolated cell with no intercell 
interference. Note that the power levels in ( fST) are less than 
what they would be in a single isolated cell. The effective 
bandwidth of each user is the same as for a single isolated 
cell, under SCP, but the power consumption is reduced by the 
factor (1 + e), which corresponds to a power gain from having 
both BSs involved in the beamforming, instead of just one. 

B. Effective interference 

The optimal SINR expressions in ( |42] l, ( |48] l, (|54]) are striking 
in how they capture the effect of interference for the three 
different beamformers. Indeed, they supply a simple "effective 
interference" characterization, which can be used to directly 
check if a particular target SINR can be achieved. 

It is natural to try and compare the schemes directly using 
the limiting SINR expressions. This is accomplished in the 
following theorem, where SNR denotes 

Theorem 4. Let "tscp^"fcBfT^MCP denote the SINRs under 
SCP, CBf and MCP, respectively. Then 



Iscp < IcBf < Imcp- 



(56) 



At signal to noise ratio SNR and interference level e, denote 
the effective interference at target SINR 7 by 



I,ff{SNR,e,-f) 



in the scenarios where the ZF beamfomier exists. 



/3| 




SNR 
1+7 


h eSNR^ 


SCP 






SNR 
1+7 


1 <iSNR\ 
1+^7 I 


CBf 


/3| 




SNR 
1+7 


1 eSNR] 
1+7 1 


MCP 
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Then the feasibility of ^ in the case of SCP, or CBf, is equiva- 
lent to satisfaction of the inequality - — > 7, and 



in the MCP case, it is equivalent to 



Ieff{SNR,e,^) 
{l + e)SNR 



> 7- 



%(5iVi?,e,7) 

Proof: Follows closely that of Proposition 3.2 in lilOi . ■ 
We note here the close parallel with the effective interfer- 
ence arising in the large system analysis of linear UL multiuser 
receivers ifTOl . This is due to the underlying UL-DL duality. 

C. Asymptotically optimal cell loading 

The above theorems characterize the optimal SINR for 
fixed cell loading /3 under SCP, CBf, and MCP, respectively. 
Our next step is to characterize this optimum loading: this 
determines the optimal number of users to serve given a 
number of antennas at the BS. This corresponds to finding the 
/? that maximizes the normalized (by the number of antennas) 
rate per cell r, i.e. the optimizing /3 that solves the following 
problem: 



maximize 



r = /31og(l+7*) 



(57) 



with 7* characterized by the appropriate fixed point equation 
(cf. Eqs (|43]), ^ and ([SSj). 

Proposition 1 (Characterization of the optimum /? for SCP). 
If 

e + p > 1 (58) 

then r(/3) is an increasing function. Otherwise, the optimum 
occurs at a finite fi* which may be found by a line search. 

Proof Refer to Appendix [G] ■ 

Proposition 2 (Optimal cell loading for CBf). If + t - 

2e^ — 1 > then r{j3) is an increasing function. Otherwise, 
there is a finite value of (3 at which r is maximized. 

Proof: Refer to Appendix |H] ■ 

Proposition 3 (Characterization of the optimum [3 for MCP). 



P 



(59) 



then r{/3) is an increasing function. Otherwise, the optimum 
occurs at a finite (3* which may be found by a line search. 

Proof: Comparing ( |55| l and ( |43| , we see that the former 
is the same as the latter with e + 7^ 4— p(i+t) • Performing 
this substitution in ( |58] l yields the result. ■ 
The above results define for each scheme a noise-limited 
region, in which cell loading can be increased indefinitely; 
however, this leads to ever decreasing rates per user, not to 
mention that more user channels would have to be learned. 

VI. Performance Results 

How do these schemes compare with each other and with 
other approaches from the literature? Clearly, CBf requires 
more CSI than SCP, and MCP involves much more BS 



cooperation, so it is not surprising that the SINRs are ordered 
as in (j56|. In this section, we obtain numerical results to 
provide a quantitative comparison between these schemes in 
different scenarios. Throughout this paper, we have assumed 
full re-use of time and spectrum across cells; however, intercell 
interference can be avoided altogether by applying the classic 
principle of re-use partitioning: we thus also consider in our 
simulations the time division (TD) scheme in which each 
BS is given a separate time-slot, which we shall also call 
"1/2-reuse". We also consider two forms of pure ZF in 
the context of single cell processing. SCP-ZF zero forces 
the same-cell interference, with the BS oblivious to other- 
cell interference. Generalized zero-forcing (GZF) is when the 
BSs independently zero force the interference in the two-cell 
system. Finally, in the MCP setting, the two BSs can jointly 
zero-force all the interference in the two-cell system, and we 
denote this case by "MCP-ZF". 

A. When is half re-use SCP better than CBf? 

Let Ptd be the cell loading in the half reuse scheme in 
which each BS transmits half the time. To compare with 
coordinated beamforming (a full re-use scheme), let /? := 
PcBf '■= Ptd/'^- Then the rate for the TD scheme is 



rTD(/?) = /3 log (l+7Ti? (/?)): 



(60) 



with -iToiP) 



2/3 



2P ^ 1+7J.b(/3) 



(61) 



The rate using coordinated beamforming is given by 

rcB/(/3)=/31og(l + 72;s/(/3)) (62) 



with 



ICBf 



1 



1 



1 



(63) 



Thus, rcBfiP) > TToiP) if e < 1, and rcBfiP) < rToifi) if 
e > 1. It follows that coordinated beamforming is only useful 
when e < 1; otherwise, it is better to partition the cells with a 
reuse factor of 1/2. Of course, if BS association is performed 
properly, e should be less than 1. 

B. Numerical results 



Figures |2]|4] compare the different schemes by varying the 
cell loading (3. We notice that when e is small, CBf gains little 
over SCP, but offers significant gains compared to pure ZF or 
1/2-reuse. When e is small, SCP-ZF is superior to GZF, as 
expected. When e is large, but < 1 (e.g. e = 0.8), CBf gains 
significantly over SCP, but does not gain much over 1/2-reuse, 
or to GZF (when the loading is perfectly optimized). Note that 
in this case, the relevant comparison is with GZF. 

When e < 1, CBf is always better, for appropriately selected 
13, than SCP SCP-ZF, and GZF If one were to insist on 
using pure ZF, and could choose the appropriate ZF scheme, 
and the exact optimal loading for that scheme, then it can 
get reasonably close to the performance of the optimized 
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CBf. Similarly, if one could select between SCP or 1/2- 
reuse, the performance can be quite close to that of CBf. 
The advantage of the latter is that it is universally good, if 
no joint transmission is allowed, across all levels of inter-cell 
interference. Compared to the ZF schemes, it performs better 
across a wider range of cell loadings. In large networks, it 
avoids the intractable frequency planning problem associated 
with fractional re-use schemes. When e > 1, MCP offers very 
significant gains over the 1 /2-reuse scheme (not depicted). 

In the two-cell model, MCP offers the most gain when e 
is large. Even when e is small, as in Figure |2] the gains over 
CBf, 1 /2-reuse, and single cell ZF schemes, respectively, are 
significant, and in Figure |4] they are higher still, because e is 
larger in that case. Unlike the other schemes, MCP improves 
with increasing e, but it requires significant cooperation be- 
tween BSs, including full data sharing, whose cost in terms 
of backhaul capacity is not accounted for here. 

Finally, we investigate the applicability of the asymptotic 
results to a finite system. In a first step, foT K = 3, N — 4, 
^ = 10 and e taking values in [.01; .1; .5; .8; 1], we solve the 



Fig. 3. Effect of cell loading on rate achieved for SNR = lOdB, e - 



optimization problems described in Section III for different 
independent samples of the channel and obtain the correspond- 
ing average rates. Even for such a small number of antennas, 
the large system analysis (LSA) results provide quite a good 
approximation. The results are shown in Figure |5] 

are time-consuming, partic- 



The optimizations in Section III 



ularly for the SCP case, which requires iterations between the 
optimization at the two transmitters until convergence. Thus, 
we consider applying the asymptotically optimal beamforming 
vectors from Section |V] (slightly modified so as not to break 
the per transmitter power constraint for the MCP scheme) in 
the finite system case. The results are shown in Figure [6] 

VII. Conclusions 

This paper has provided an asymptotic analysis of a two 
cell interfering network in which the number of antennas at 
the BSs, and that of users per cell, both grow large together 
Schemes that balance rates across users in the system were 
compared for three levels of cooperation, namely single cell 
processing, coordinated beamforming and multicell process- 
ing. MCP offers significant rate gains, if we can accommodate 
the additional coordination and communication between BSs, 




Fig. 4. Effect of cell loading on rate achieved for SNR = lOdB, e = .8 



not accounted for here. We characterized and compared the 
limiting SINRs of the three schemes, in particular, using the 
notion of effective interference, which can also be used to 
establish if a given SINR target is feasible. The validity of the 
obtained results was verified via Monte Carlo simulations. 

Note that we have assumed the users' channels are selected 
randomly. It is important to emphasize that our conclusions do 
not hold if users have been selected based on their channels, 
as the conditional distributions change drastically. Indeed, it is 
well known that for a large enough pool of users, with careful 
scheduling, the performance of ZF can be almost optimal. 

Appendix A 
Useful Theorems 

We start by reproducing a few lemmas, which play an 
important role in our derivations. 



Lemma 1 (Lemma 6.3.3 in 1|52J). Let p>0, AandB NxN 
matrices with B Hermitian, r € M, and q € C^. Then 

lAII 



tr 



(B- 



B + qq^ + pi 



< 



P 



(64) 



We will also be making use of results from 11531 (some 
themselves reproduced from ll54l ). Il53l derives a central limit 
theorem for the SINR at the receiver in a multiple-access 
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0.2 0.3 



Fig. 5. Average normalized cell rates for the optimization problems for 
K = 3,N = 4 and LSA rates for B = .75, both at SNR = 10 dB. 




SNR (dB) 

Fig. 6. Average normalized cell rates from using the asymptotically optimal 
beamformers to the finite system case for K = 3, N = 4 and LSA results 
for /3 = .75, both with e = 0.5. 



MIMO system, where the dimensions of the system (number 
of users and the size of the received random vector) grows 
large: the considered asymptotic regime satisfies 

K — > cx),limini — > 0, limsup — < oo, (65) 
N N 

and the quantity studied, i.e. the UL SINR, is equal to 



r,^^y" (YY^ 



(66) 



where the sequence of matrices S(i^) 
given by 



yiK)YiK) 



IS 



\ / n=l,k=a 



N,K 



-Wnk 



= l,k=0 

(67) 



and is such that the following assumptions hold: 
Al: the complex r.v.'s {wnk n > 1, k > 0) are i.i.d., Ewio = 
0, Ewlo = 0, E|wioP = 1 and E|wio|^ < ocj^ 
A2: the variance profile is such that there exists a real number 
o'max < OO such that 



sup max \S-nk(K)\ < d-„ 

X>1 l<n<Nfi<k<k 



(68) 



*For the channel model considered here, the entries are i.i.d. CM{0, 1), so 
they clearly satisfy this condition. 



The dual UL SINR expressions in the present paper for Afej's 
in each cell held at some constant value, will be of the 
same form as (|66|, and the matrices considered will satisfy 
assumptions Al and A2. In fact, they constitute a special case 
of the above model, since we consider the asymptotic regime 
where ^ = c as N ^ oo, and the variance profiles of 
the random matrices considered (see ( |67] i) are either scaled 
all ones matrices, or obtained by the regular sampling of 
a piece-wise continuous function. Note that the expressions 
for asymptotic SINR on the dual UL's may equivalently be 
obtained from earlier results in the literature, 1 11 1 for example. 

We thus reproduce below the most relevant results from 
1531. 



Theorem 5 (Parts 1 and 3 of Theorem 1 in [53]). The 

following statements hold true. 

• Let ^o'^fe(^) :l<n<iV, l<fc< be a sequence 
of arrays of real numbers and consider matrices T>k{K) 



Bk{k) = diag {~aUk),..., dlAk) ),0<k<K 



The system of N functional equations 

1 



n,K 



{■A 



for I < n < N and where 



k=l l+^,rD,{K)T^(z)' 



(69) 



(70) 



(71) 



admits a unique solution T among the diagonal matrices 
for which the ^ belong to class iSj^ Moreover, the 
functions admit analytical continuations over C — [0, oo) 
which are real and positive for z € (— oo, 0). 
Assume Al and A2 hold true. Consider the sequence of 
random matrices Y^ {K)y (K)^ , where Y„fe = ^j=Wnk- 

For every sequence Sk of N x N diagonal matrices with 



T^(z) =diag[t^^j^{z),...,tff j^{z) 



sup I 

K 



'A' I 



< OO 



(72) 



the following limits hold true almost surely (a.s.): 



lim 



K- 



K 



K 



(Q^(z)-T,^(^)) = 0,Vze 



(73) 



where Q^(z) denotes the resolvent of {K)Y {K)^ , i.e. 
the N X N matrix defined by 

Q^(z) = (Y{k)Y{k)" - zl^) . (74) 



_ Corollary 4. Assume Al and A2 hold true. Let = 



j^tr(Sj,Cl\{-p)) where 



'K' 



— ifrS^^ (T^(— p)), and be as given by Theorem^ 
Then 



a.s. 



(75) 



'A complex function t{z) belongs to class S if t(z) is analytical in the 
upper half plane C+ = {im(z) > 0}, if t{z) S C+ for all z S C+, and if 



im2: \t(z)\ is bounded over the upper half plane C 



+ ■ 
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Proof: Differentiating ( |73] l with respect to z, we get where 



Jim -^tr(S^Q|(z)) - -trS^^ (T^(^)) = 



K-s-oo K 



K 



Vz e C- 



(76) 



^ (T^(z)) exists since j^{z) admit analytical continua- 
tions over the range considered. This yields the result. ■ 



Theorem 6 (Theorem 2 in 1531). Lef = 

-ifr ^Do(-ftr)T^(— p)^ where p G M_|-, and w ^/ve« fey 



Theorem |5] AiiMme A7 ant/ A2 then 



r 



r 



K 



a.i. 



(77) 



Proposition 4 ( 1531 ). Introduce the r.v.'s Ui — -irrDoQD;Q, 

/or Q < I < K. The Ui satisfy the following system of 
equations: 



^ 1 
Ui = Y^ cikUk + -^tr-DoDiT^ + Q, < / < X, 



where 



Clk 



1 prBi-DkTi-pf 



K 



(l + ifrDfcT(-p) 



(78) 



(79) 



and the perturbations ti satisfy E |ei| < CK ^1"^, where C is 
independent of I. 

Appendix B 
Useful Asymptotic expressions 

In this appendix, we derive asymptotic expressions for 
quantities of interest in the large asymptotic analysis of SCE, 
CBf and MCP for the considered two cell symmetric channel 
model. This relies on the application of the theorems in 
Appendix |A 
in Theorem 



[noting that for our model the T diagonal matrix 
5] collapses to one or two variables, as many of 



its entries will be equal. 



A. SCP 



Clearly, our Rayleigh channel model satisfies the conditions 
in Theorem [6] so that for j's fixed for all users in cell j at 
boundecrj A j , the quantity of interest for any user is 



^'rh,,,,, ( ^ h0,,,,h.',,., + I I h^^,,. (80) 



N 



k'^k 



From Theorem |6] this quantity will converge a.s. as N,K 



oo with 



scp,UL _ Ajtscp(-1, Aj), 



(81) 



'"This is required for the results in Appendix |a| to be applicable. On the 
other hand, if this was not the case, in the studiedoptimization problem, the 
dual objective would be unbounded and the primal unfeasible. 

"The T matrix in the theorem is simply a scaled identity with scaling 
factor converging to tsQp(~l,\j), as given below. 



tscp{z-,\j) 



Thus, rySCP,UL j-;^^.) satisfies 



7 



SCP,UL 



A, 



1 



(82) 



(83) 



Applying ( |73l ), we can show that 



1 

— t 

N 



tscp{-p, Aj). 



Furthermore, applying Corollary [4j we get that 

\ -2 



1 

TV 



tr 



A 



(84) 



(85) 



converges a.s. as N,K ^ oo, ^ — > /3, to 
d 



dp 



tscpi-p, Aj) 



p=i 



(86) 



From ([82]l, tscp{-p, Aj) = — — , so that 



dp^scp{~P,Xj) = - 



{i+\jtscp(-p,\jyf 



(87) 



B. CBf 

Here too, we may apply Theorem |6] when A^ j 's for all users 
in cell i fixed at Xj, and the values of /ij also held constant in 
both cells. For any user in cell j, we will need to characterize 



N 



(88) 



where 



H 



^h|(,,,hfe,,,, + |5:h^_._^h,,j,. (89) 

k' 



N 



k'=tk 



From Theorem |6] this quantity will converge a.s. as iV, iiT — > 
oo with f ^ ^, tcji2] 

^(^BL.UL (^^., A,, Aj) = A,<cs/(-/ij, A,, Aj), (90) 

where 

tcB/(z, Aj, Aj) 

— . (91) 



-z + 



(l+AjtcB/(z,Aj,A3)) (l+eAjtcB/(2,Aj,Aj)) 



This uses the fact that 



1 and ^ — > /9 as ii" — > oo. 
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Thus, -yCBf-.UL (^^.^x.j,Xj) satisfies 



n. + + 



Applying ( [73] ), we can show that 

tcBf{-P: Aj , Aj). 
Moreover, by Corollary |4] we can also show that 

\j' = l k' I 



(93) 



(94) 



converges a.s. as iV, if — > oo, ^ — >■ /3, to 



dp 



tcBf{-p, Aj, Aj) 



(95) 



From (|9T) 



dp 



tcBfi^P, Aj , Aj) 



-tcBf{-p, Aj , Aj) 



(l+AjtcB/(-p,A,,A5)) (l + eAjtcB/(-p,A,,A3)) 



(96) 



C. MCP 



For Afc j's in cell j equal to a constant Xj and the pj's 
also fixed held fixed, for any user in cell j, we will need an 
asymptotic expression for (cf. Eq. (|34|l) 



k' 



N 



k'^k 



life J- • 



(97) 



Introducing the short-hand = (pj, pj, Aj, Aj) , applying 
Theorem |6] we can show that ( |97] l converges a.s. as N,K 
oo with /3, to ry^'^^P-^^ (^^ij,fij,Xj) equal to 



7i (»?) = Aj ' 



Pj 



M.7 



(98) 



where ti^MCp{z-,r]) and t2,MCp[z,r]) are given by ( |99| ). 

For a user in cell j, one can verify that ( |97] i (replace j by 
j and vice versa) converges to 

MCP,UL , N , / il.MCp(-l, ''7) , ^2,MCP(-1,''7) 

72 (»7) = Aj e- 



Pj 



pj 



(100) 



Plugging in ( |98| l and ( |100[ ) in the expressions for 
ii,j\/cp(-l,»7) and t2,MCp(-l, »7), (cf. Q), we get 

1 

ii,j\/cp(-l,i7) = 



1 



^7 



MCPML I 7 1 T"; MCP,UL f 7 
l+7i iv) 1+72 (^) 



(92) i2,MCp(-l,T7) 



l + 7l (»?) 1+72 (^) 



(101) 



Now using ( fToTj l in (|98) and (fTOO|, 7i ' (ry) and 

MCP,UL I X ^Ti f 

72 ^ill satisfy 

MCP,UL I N 

7i {m 

= A 



;3A, ^g^j 



/3eA, 



l+7i (v) 1+72 ■ ('?) 



MCP,UL 
72 



= A5 



Pj -T ^^^MCPML 

iv) 



13^-, 



fix, ^g^j 



/9eA, 



Pj ^ YWUTTWI: 

Now define 



/3A. 



l+7i (V) 1+72 ('?) 



Dfci 



Ml 



fJ-2 



Onxn e 



Dfc,2 = 

Applying ( |73] l, we can show that 



A 



(102) 

(103) 
(104) 



a.s MCP.UL I \ 

— > Ij iv) ■ 

We also need to characterize 

2 



1 

TV' 



-trD 



j'=i k' 



(105) 



(106) 



for the special case where p^ = pj = p and Aj = Aj = A, 
i.e. T] = rj^y,^ = [p, p, A, A]. In this case, 

h.AICPi-P^Vsym.) = h.MC p{- Pi V syrn) = tMCPi-PiVsym) 

^ —■ (107) 



Thus, 
d 



(l+e)/3A 

tJ, + (l + e)XtMCP(.-P,Vsym) 



, . iMCpi-PiTlgy^) 
^tMCP(-p,Vsym) (l+e)/3AM 

[^l + {l + e)XtMCP(-p,Vsym)Y 



(108) 
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KMCpiz,ri) - — , 

-z + + 

l + 7r-ti,MCf>(z,'7) + 7fT^3*2,jv/Cp(z,»7) l+« 77^ ti,af CP (z,»?) + T7^ A -ta.MCpCz.r?) 

t2.MCp{z, rj) = ^ . (99) 

-z H h 

l+7f-^jtl,MCp(z,I7)+<;^t2,MCp(2,»?) l+fT>^-jti_^MCp(z,n)+jr^t2.MCp(z,n) 



and, by Corollary |4] ( |106[ l converges a.s. as A^, if — > cjo with 

f ^/3to 



(l + e)A d 



fi dp 

Finally, we will need to characterize 
1 



sy7n ) 



N 



trDfe jAD 



(109) 



(110) 



with A as defined in ( |194| i, where in Ai and A2 in 1, 
Dfe/ 2 are equal to A for all k < K. For any sequence of 
diagonal matrices S with bounded diagonal entries, and taking 
into account the fact that the D/j.^'s are equal for all users in 
the same cell, define Vj (S) = ^^trSADijA, for j ~ 1,2. 
Applying Proposition |4j we get that 



Vi (S) = ^ 



/3- 



{l + {\ + e)\tMCp) 
{l + {l + e)\tMCp)' 



MCP 



N 



V2 (S) = 13- 



trSDi, 

'^'^'-MCP 



A2 



'l + {l + e)\tMCp) 

1 + (1 + ()\tMCp) 



iVi (S) 
2 ^2 (S) 

^^1 (s) 

2 ^2 (S) 



(111) 



t 



^ti-SDi . 

TV 



(112) 



where 771 and 772 satisfy E|77j| < CN^^/"^, j — 1,2, and 
tMCP ~ tMCp{^, 1, for some constant C. 

For S = Di.i, we obtain ( |1 13[ ) and ( |1 14| i, where and (/)2 
denote deviation terms such that E\(j}j\ < CN^^/"^. 

Note that by definition, Vi (Di,2) = V2 (Di^i), and 

1^2(Di,2) = Fi(Di,i). 

Appendix C 
a simple monotonicity result 

Let /(A) be a standard interference function for the UL, in 
the sense of Yates ||27l, where we denote the transmit power 
vector by A. Suppose it is of the form 



h{\)^lM\) k^l,2,...,K 



(115) 



where is the SINR target for user k, and F{X) is a standard 
interference function. The vector A is called /efli/Z7Ze if 



Xk>lkFk{X) k = l,2,...,K. 



(116) 



It is shown in fZTl, Theorem 1, that if a feasible solution 
exists, then function / has a unique, positive, fixed point. 1,271 
also shows that, starting at any power vector A, the iterative 
power control /"(A), n — 1,2,... converges to it. Using the 
fixed-point powers, all users achieve exactly their SINR target. 

Now consider two different vectors of SINR targets 7'^' 
and 7^^-', and let A'^-* and A'^-* denote the corresponding fixed 
points. Lemma |2] is a simple corollary of Lemma 1 in ||27I . 

Lemma 2. 7/7(1) < 7(2) then a'^^ < A^^'. 

Proof: A'2) is feasible for the power control problem with 
SINR targets given by 7(1). By 1271, Lemma 1, /" (a^^^), 
n = 1,2, .. . is a monotone decreasing sequence of feasible 
power vectors that converges to A'-^'. ■ 

Appendix D 
Proof of Theorem[T] 

Throughout this section, let 



A 



1=1 



where A > will be defined later 



(117) 



(118) 



(119) 



Asymptotic analysis of the dual problem 

We start by considering the dual problem at each of the 
BSs: this will yield the asymptotically optimal dual variables 
and beamforming directions in both dual and primal problems. 
Note that even though the value of the dual objective function 
in one cell depends on the primal beamforming decisions 
through the (Tfc j S in ( [T4] i, the optimal dual variables themselves 
are fully determined by the constraints and will be the unique 
strictly positive solutions to ([TSj. Since the analysis is identical 
in both cells, without loss of generality, assume the cell index 

j is j = 1. Assume also that j3 



7 



1+7 



£7 



< 1. As noted 



in Section [V] we cannot immediately apply standard large 
system analysis to ([TtJ as optimal dual variables A^.i's are 
not independent of the channel vectors (hfc_i,i)^^j^. 

Rather than directly analyze the asymptotic performance of 
the SCP system using optimal Afe.i's, consider any constant 
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V, (Di.i) = tljcp>? (1 + eY 



V2 (Di,i) = tl,cp>? (1 + eY 



^(l+(l+e)AtMCJ.)' 



1-/3 



1-/3 



(l+e)2A2tJ 



(l + (l + £)At„, 



1-/3 



^MCP 

{l+il+e)\tMCpY 



It 

(l+e)2 



(l + (l + £)AtMCp) 



(l+(l+e)AtMCp)'' 



(113) 



(114) 



A > (later, we will assign it a particular value, but for now 
it is arbitrary) and consider the large system regime in which 
all users have the same transmit UL power of A/iV, as iV, iiT ^■ 
00, with fixed ratio KjN = /3. This is a virtual UL with noise 
of unit power, so the SINR for user fc is given by 



\_. 



ife,i,i 



I 



-E 

k'^k 



hfc^,i,ihfc'a,i 



1^ 

ifc.i.i- 



(120) 



Appendix |B-A| shows how this converges a.s. to a constant in 

7 + 5 



the considered asymptotic regime. 

In particular, if S is small enough so that /3 



1 + 7 + (5 

and we assign the following particular value to A: 



< 1, 



1^(5 J 



+7+S 



(121) 



which we denote by X{S) to make explicit its dependence on 
the parameter 6, then ( |120[ ) will converge a.s. to j+6. Theorem 
3 in II53I shows that for this (suboptimal) system, the value of 
(fnok is 



Vfc = l,2,...,if, (122) 




where jk{6) denotes the dual UL SINR of user k under this 
suboptimal power allocation, and the last term on the right 
hand side is 1 /^/N times a r.v. that converges weakly to a zero- 
mean Gaussian distributed r.v. This holds since ^ ^ /3 < 00, 
for the considered channel model and in the notation of the 
theorem, and 0^ will converge a.s. to bounded limits. 

Denote by X{S) the vector of UL powers in this suboptimal 
system, where all entries have the same value X{S). For (5 > 0, 
it follows from ( |122| l that for any particular k, P(7/c(— (5) > 
7 or 7fc((5) < 7) decays to exponentially in N (it is a large 
deviation event). Applying the union bound, we obtain that 



7fc(-<5) < 7 < 7feW Vfc = l,2,...i^ 



(123) 



will hold with probability tending to 1 as TV t 00. Let 7(— (5) 
and 7 ((5) denote the vectors grouping the left and right-hand 
sides of ( |123| l, respectively. 

Denote by A the vector of optimal UL powers for the dual 
problem, which achieves SINR of 7 for each user On the 
other hand, A(— 5) and X{6) (defined above) are the vectors 
of optimal UL powers for the virtual UL problem with target 
SINR vectors 7 (—(5) and 7 ((5), respectively (instead of the 



all-7 vector in the original dual problem). Thus, by Lemma [2] 
in Appendix |C] 



A(-(5) < A < X{6) 



(124) 



will hold with probability tending to 1 as iV f co- Since 
this holds for any sufficiently small 6 > 0, the empirical 
distribution of the components of the optimal A will converge 
weakly to the constant A, given in ([39|. Let ([39| provide the 
particular value of A in the rest of this section. This establishes 
the asymptotic optimality of having the Afcj 's in both cells all 
equal to A, and consequently that of the beamforming vectors 



= Afc.,h 



H 



kJ'^^k,j,j- 



(125) 



By UL-DL duality and from the KKT conditions, these are, 
up to a scale factor, also the optimal DL beamforming vectors 
|29l. Thus, this analysis shows that, asymptotically, the DL 
beamforming directions in one cell do not depend on the 
beamforming directions used in the other cell, although it 
is clear that the optimal DL power levels (which modulate 
the beamforming directions) will depend on the power levels 
used in the other cell, even in the limit. Finally, also note 
that the dual objective function value SfcLi '^k j ^'^ 
upper bound on the primal objective function, i.e. on the total 
transmit power in the cell; denote the latter by Pj for cell j. 

Asymptotic analysis of the primal problems 

We now turn to the DL primal problems, and fix the 
beamforming directions in both cells to be those given by 
( |125| l. Thus, only the DL power levels pkj in ( [T9] l need to be 
determined. These must satisfy, for all users in both cells. 



Pk'j 

k'^k N 
1 |ht,j,jv 



"k'j I 



l2 
'kj\ 



K 
k' = l 



h 



kJJ^k'j 



N 



w^k'jr 



(126) 



(127) 



Assuming feasibility, one can verify using similar standard 
interference function arguments as for the dual problem, that 
the set of equations in both cells will have a unique power 
minimizing solution |[26l . 

As was the case with asymptotic analysis of the dual 
problem, it is easier to first fix the DL powers to constants 
and study the resulting limiting regime. Thus, assume pk/s in 
cell j are held fixed at a constant value pj j = 1,2, and note 
that this implies that Pj = jSpj. After analyzing this limiting 
regime, we optimize the choice of the constants, pi and p2 
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and finally show that the optimal constants are asymptotically 
optimal with respect to the primal optimization problem (j6|. 

Under the regime in which pkj — Pj for all k in cell j, the 
following lemmas hold. 

Lemma 3. With -w^j — Akjh^jj, such that A = i 3 i ' 



the following holds, as K, N 00, ^ — (3. 



max 

j=l,2,k<K 



max 

j=l,2,k,l<K 



(i+7r 
2 



w 



max 

j=l,2,k,l<KA^k 



|h/c,jjwy| 1 



> (128) 

> (129) 
0. (130) 



Proof: Applying Lemma 5.1 in ll55l . 



max 

j=l,2,k<K 



N 



1 

N~ 



1 2 - -trDfe,,,,A^,^ 



max 

j=l,2,k<K 



^0 (131) 
' 0, (132) 

wher^*^ Dfcj.j — ^I- Later on in this section, 'Dk.jJ = ^^I- 
Now applying Lemma [T] twice to ( |131| l and once in ( |132| i, we 
get 



max ^llwfejlp - ^trDfe.jjA^ 

j=l,2,k<K N N ■' 

max — hi, i iWki trDi, „ , A, 

j=i.,2.k<K N ^ N ^ 



>0 (133) 
0. (134) 



Now consider the interference terms, with i — 1,2, and 



A 



(135) 



Two different cases arise, depending on whether i = j or 
i = j. In the latter case, we may apply Lemma 5.1 in ll55l to 
the numerator of the right-hand side of ( |135| l. 



)2 



0. (136) 



Applying Lemma 5.1 in f55l once again yields. 



max 

k,l<K 



^l3J-^lJ^k,j.j-^lji^ 



H 



N 



-trDiJjA;jDfc_^jA;j 



Finally applying Lemma [T] twice, we get 



^"■DiJjAjjDfc -A,j - ^trD, , ^ A, , ^ A, 



N 



N 



1 



< — D, -A, J) 



< 



N 
N ' 



< 



I^'JjII II A' j II ll^fcjJ 



iV 



trDiJjA/jDfcjjAj - -trDfejjAjD^jjAj 



1 

7 



(138) 



(139) 



Thus, 



max 

k,l<K 



0. (140) 



When i in ( |135| l is the same as j, we cannot apply Lemma 
5.1 in 155 1 directly, since A; ^ and h^jj are not independent. 
Thus, we apply the matrix inversion lemma first to get 



(141) 



(1 + w^k,j,jAk,i,jh^.jj 
We can now consider the numerator in ( |141| i and show that 



max 

k,l<K,(k^l) 



N 



hfc J J- Afc ^ ; J j h; J- J A; h|^^- J 



N 



0. (142) 



We can also show that 
A 



max 

k,l<K,{k,yil) 



N 



N 



:tTT>k,j,jAj 



^ 0. 
(143) 



The proof of the lemma is concluded by us ing the limits 



B-A 



with A, = A, 



of the trace terms as derived in Appendix 
noting that for this specific case, all the 'Dk.j,j and Dj^. ^ j 
matrices are equal, and are simply scaled identities, so that 



-^trT)k,j.jAj = — trAj 7 



1 

N 
1 

N 



A 

trA, 

A 



(144) 



-trA 



2 a.s. 



BX 



Ai_ 



— trDfc , ,-A,D, , ,-A, = —trA? ^ ^ 



1 



0. (137) -trDfej,j.A,Di,,,,A, 



N"' ' 
eA2 



iV 



N 



-trA 



2 Q.S; 



(1 + 7)" 

2 



7 



(145) 
(146) 

(147) 



(1+7)" 



'^We introduce the Dj; j ^ matrices to allow for more general formulations. 
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Lemma 4. With \Vkj = Akjhf^jj, such that X = i ^ , 
and with pkj = pj for k = \, . . . ,K, j — 1,2, it follows that, 
with probability 1, 

al^^a^ + ePj, (148) 
Pj denotes the total transmit power of BS j, i.e. Pj = l3pj. 
Proof: 

Since DL pkjs are all fixed to the same value, pj, we have 



2 „2 .Pj \^ |hfcjjW;,j|2 

(Tl. ; = (T H > 



N 



(149) 



k=i " ''=^■11 

Applying ( fT29] l of Lemma [3j for fc = 1, . . . , if, j = 1, 2, 



a^ + ePj + o{l). 



(150) 



We conclude from Lemma H that cr^ j is asymptotically 
independent of the user index k, in the regime considered 
(Pkj = Pj) and therefore write the limiting value as ct|. This 
deals with the asymptotics of the RHS of ( |127| l, and we now 
provide two lemmas to deal with those of the RHS of ( |126| l. 

Lemma 5. With wt , = At oh? ■ -, such that A — — J -v , 

■I ■■■> i^jj ^~Py^ 

and with p^j — pj for k — 1, . . . , K, j = 1,2. The following 
holds for any k' in cell j: 



N 

k'^k 



llWfc'jII^ 



- (3 



(1 + 7)' 



(151) 



as K,N ^ oo, ^ 
Proof: 



/3; A- = I3P,. 



^ N 

k'^k 

K - 1 



Pj ^ \tlkJ,jWk'j 



N 



k'^k 



N (1+7)2 



0(1), 



(152) 



where we used ( |130| l in Lemma [3] Noting that as i^, — > oo, 



^^y^ — ?► 13 completes the proof. 



To find the minimal pair of constants {pi,p2) for the two 
cells, we therefore solve the following set of equations. 



Pi 



P2 



7 



(1+7)^ 

7 



_M1 



ePpi 



Pi 



P2 



to obtain 



(1+7)^ L 



Pl^P2=P 



(1 + 7)' 
(1 + 7)' 



1 - 



(1+7) 



e/37 



This implies that 



Pi=P2 = P:= pp. 



(153) 



(154) 



(155) 



(156) 



Note that such a choice of transmit powers is guaranteed to 



meet the SINR constraints as K, N 



K 

N 



We confirm the asymptotic optimality of this deterministic 
power allocation, together with the beamforming directions 



found from analysis of the dual problem, by verifying that 
the duality gap is tending to zero in both cells. Indeed, in 
cell j, the primal objective value converges to P, whereas the 
dual objective value converges to /3A {a^ + eP). Recalling the 
definition of A, one can easily verify that these two quantities 
are the same, thereby completing the proof. 

Equation \155\ shows that the coupled primal prob- 
lems have a solution when effective bandwidth condition 
/3 (yY^ +£7) < 1 is satisfied. Conversely, if this condition 
does not hold, there is asymptotically no feasible solution to 
the coupled primal problems, in the limit as A t 00. The latter 
observation follows from monotonicity: The optimal powers 
are increasing functions of 7, but as the denominator of \155\ 
decreases to zero, the optimal total power from either BS, P, 
must diverge to infinity, and for higher values of 7 there can 
be no feasible solution. 

Appendix E 
Proof of Theorem|2] 

Throughout this section, let 



A 



2 K 



A.= (l+^EEh5> 



l-i,j 



(157) 



i=l 1 = 1 
2 



A^.^=|I+]^E E hf^Mhu,,) (158) 



2 



Afe,,,fc',,',,= |i+^E E K,j^^ 



{t.l)^(j,k),{j',k') 



(159) 



where A > will be defined later 

Moreover, the following lemma will be useful in proving 
the theorem. 

Lemma 6. Let {^1,^2) satisfy < fii < 2, i — 1,2, and 
define the function F{Xi,X2) by 

F(Ai,A2) = (Fi(Ai,A2),F2(Ai,A2)), 

^;(Ai, A2) = 7 + A, + Xj j = l, 2. 

(160) 



Then 

i) F is an interference function 1127V . 

/ 7 ^7 \ 

ii) If B 1 < 1 then there exists a unique 

y 1 + 7 1 + £7 / 
solution to the fixed point equation 

CXi,X2) = FCXi,~X2)- (161) 

Proof: (i) is easily verified. For (ii), let the value A be 

27 



A = 



(162) 



i-/?(t^^ + t^)' 

Since ^i,/i2 < 2, it is easy to verify that 

F{^, A) < (A, A) componentwise (163) 
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and hence the inequality (Ai,A2) > F(Ai, A2) has a feasible 
solution. It follows from Theorem 1 in |27| that ( |160| l has a 
unique fixed point (Ai, A2). The latter depends on {fii, fi2). ■ 
Asymptotic analysis of the dual problem 



7 



£7 



be 



Assume that /? 1 < 1, let (^i,^i2 

Vl+J l + e7/ 
feasible for the dual (|22J i.e. /ii,/i2 > and /ii + //2 = 2. 

Consider the suboptimal dual power vector that assigns power 



level Xj/N to all the users in cell j. 
The dual UL SINR is equal to 



j = 1,2, respectively. 



For the purpose of obtaining a contradiction, let jli be such 
that Fi {p.1 -5,fLi+S)>0 for all 5 > 0, and let /i2 = 2 - /2i. 
We will assume < /2i < 2{^ Roughly speaking, there is 
non-negligible probability that ni will be close to /2i when N 
is large, along the subsequence. Indeed, for 5 > 0, let us define 
Bi{d) be the event that fii e {fli ~ 6,fli + S), then by the 
second Borel-Cantelli lemma, event Bi{S) will occur infinitely 
often. Due to the feasibility constraint that /ii + /i2 = 2, we 
can equivalently write 



N 



2 

E 



N 



(164) 



We show in Appendix |B -B | that this quantity converges a.s. to 
a constant equal to AjtcB/(— Mji %) (cf- Eq. (|9T|). 

Letting (Ai,A2) be the unique solution to ( |161| l, 
^j^CB/(— A'j, Aj) will be equal to 7 for j — 1,2. Thus, 
with this set of suboptimal dual UL powers, all SINR's 
(cf. Eq. ( |164[ )) converge to 7. Note that the condition that 
< A*; < 2, i = 1,2, is weaker than the condition that 
{ill, ^2) satisfy the dual feasibility condition that /ii +/i2 = 2, 
but it certainly includes this condition. 

We now arrive at a particular choice of (/ii, /i2) by solving 
the following optimization problem: 



/3(T^(Ai + A2) (165) 
s.t. fii + /i2 = 2 (166) 
(Ai, A2) is the unique fixed point of ( |161| l. (167) 

Adding up the two equations in ( |161| ), and taking account of 
(|166|), we obtain 



1 



/3 



1+7 
27 + ^e 



(Ai + A2) 

A1A2 ^ 
Ai + 7eA2 A2 + 7eAi 



A1A2 



(168) 



Let t — Ai + A2 for (Ai,A2) optimal. For fixed t, the 
(^11^2) that maximize the RHS of ( |168[ ) can be verified to be 
(t/2,t/2). The optimal Ai and A2 must thus be equal; we de- 
note the common value by A. ( 161| l then implies jii ^ 112 ~ 1- 
As a result, A is as given in ( |45] l. We conclude that the optimal 
choice (with respect to the optimization problem ( |165| l-( [T67| i) 
is to take {^,1,^2) = (1, 1). 

Note also that if we fix (/ii,/i2) — (1,1) and use the 
corresponding deterministic (Ai,A2) = (A, A) (with A satis- 
fying ( |45] )) to generate dual UL powers in each cell, then the 
dual objective function will converge to the solution of the 
optimization problem ( |165| l-( [T67l ). 

Now consider the optimal dual variables (/ii,/i2), and the 
vector A of {Xkj)^Li^^^i, where optimality here refers to the 
dual CBf optimization problem ( |22| i. Due to the dual feasibility 
constraints, the sequence of (/ii, 112) is contained in a compact 
set and so the probabihty distribution function of /^i, f[^\ 



that F, 



forms a tight sequence 115611 . Let Fi denote a limit point, so 



^1 — (5<^i <//i+(5, ^2 — '5<^2 <A'2+<5 (5>0 

p.l+S<p,i<p.i—S,p.2+S<fl2<P-2—S S<0. 

(169) 

• To compare the performance of the optimal scheme with 
a deterministic power scheme, we modify the SINR target in 
( |160| l from 7 to 7 + 5 for some small constant S (positive or 
negative), and we replace {p,i,p,2) by {p,i + S,fl2 + S). The 
latter change violates the constraint that the sum of the /is 
should be 2, but it nonetheless provides a valid pair of noise 
values for a virtual UL. Provided that \S\ is sufficiently small, 
the condition < jli + S < 2 i = 1,2, will be met. 

Denote the corresponding solution to ( |161| l (with (/ii,/i2) 
replaced by (^1 + S,fl2 + S)) by (Ai((5), A2((5)), and let X{S) 
denote the vector of UL powers, where cell j users use power 
level Xj{S). Note that {Xi{S),X2{S)) (and hence \{S)) depends 



on (/2i, /I2), and S. The analysis in Appendix B-B shows that 
this power allocation will asymptotically achieve a SINR of 
7 + (5 for all users, under external noise levels of /2i + S, 
and /i2 + 6, respectively, and for finite N, let it achieve SINR 
lkj{5) for user k in cell j. 

As in the SCP case. Theorem 3 in ll53l can be applied to 
show that for this (suboptimal) system, the value of (|164[) is 



lko{5) = ^ + 6 + 



where the last term on the right hand side is 1/viV times a 
r.v. that tends weakly to a zero-mean Gaussian rv. This holds 
since ^ /? < 00, for the considered channel model and in 
the notation of the theorem, and Qj^ will converge a.s. to 
bounded limits, though these will be different from the ones 
in the SCP case. As in the proof of Theorem [T] in Appendix 
[D] a union bound may be applied to show that 




Vfc-l,2,...,if,j-l,2 
(170) 



(171) 



Fi along a convergent subsequence. 



Ik, {~5) <l<lkA5) yk = 1,2,..., K,j ^1,2 

will hold whenever Bi{6) occurs, in the limit as N ^ 00. 

Now compare the dual optimal power levels with this 
suboptimal, deterministic power allocation, on the event that 
Bi{S) occurs. Consider the case that S > 0: Noting ( |169| l, 
we see that, on this event, the dual UL powers required to 
achieve jk.ji^) will only decrease if we replace the enhanced 
noise levels, {p,i + S,fl2 + S), by (/ii, 112)- But by ( |170| l and 
the monotonicity result in Lemma |2] in Appendix [C] it follows 
that these decreased UL powers must upper bound A since the 
latter is the optimal vector of powers in each cell to achieve 

'''The cases /Ii = or /ii = 2 can be considered separately, in a similar 
manner, which we mention at the conclusion of the proof. 
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the all-7 SINR vector under noise levels (/^i,/i2). Similarly, 
consider the case ^ < 0: On the event Bi{6), the dual UL 
powers required to achieve j^jiS), will only increase if we 
replace the lower noise levels, (/li + 6,jl2 + S), by {111,112), 
but by ( |170[ ) and the monotonicity result in Lemma |2] in 
Appendix |C] it follows that these increased UL powers lower 
bound A. We conclude that for N sufficiently large, and on 
the event Bi{S), we have that 



Ai(-^) < Ai < Ai(<5) 
A2(-^) < A2 < X2{S) 



(172) 
(173) 



and hence that the optimal dual objective value lies in the 
interval {Pa'^{Xi{-S) + X2i-d)), l3a^{Xi{S) + X2id))). But S 
can be taken arbitrarily small, so the left and right endpoints 
of the interval can be made arbitrarily close to /3cr^(Ai + A2), 
where (Ai, A2) are the corresponding deterministic UL powers 
when we set 6 = 0: These are the unique solution to ( |161[ ) 
for the given (/li,/i2). Thus, when N is sufficiently large, 
along the chosen subsequence, and on the event Bi{S) (which 
occurs infinitely often) we have the dual objective value getting 
as close as we like to the value f3a'^{Xi + A2), which is a 
feasible value to the optimization problem ( |165| l-( [T67| ). Thus 
if (/2i,/l2) 7^ (Ijl) ™d < /ii < 2, then we get a 
contradiction of the optimality of the optimal dual power 
levels, whenever Bi {6) occurs, since we can beat the purported 
optimal value using the deterministic UL power levels obtained 
for (/ii,/i2) — (1,1)^ i-S- ( |45] l: this is the case since as 
discussed earlier the corresponding dual objective converges to 
the solution of ( |165[ )-( fT67] i. A similar argument can be used to 
show that {fli, p,2) — (0, 2) or (/ii, /i2) — (2, 0) lead to similar 
contradictions. Since these are contradictions, it follows that 
(Mi,M2) = (Ijl) and hence all the mass of the distribution 
Fi must be concentrated at (1, 1). 

We can conclude from the above analysis that for any 
positive 5, the event Bi{5) will occur with probability tending 
to 1 as tends to infinity, and therefore that ( |172| i-( [T73| l hold 
with probability also tending to 1. Taking 5 to zero, we obtain 
that the empirical distributions of both the Ai^s and the A2A;S 
tend to the same constant, namely A, as N tends to infinity, 
where A given in ( |45] i. 

We conclude that the asymptotically optimal dual variables 
are /ii = /i2 = 1 and all A^jS converge to A as given 
in ( |45] l. Thus the optimal downlink beamforming vectors 
asymptotically point in the directions of the vectors given in 
( |47j l. It remains to find the optimal power levels, pkj, to use 
for each user. 



Lemma 7. With Wk-j = Ak.jh^jy X satisfying ( |45] l 



max 

j = l,2.k<K 



max 

j = l,2,k,l<K 



1-/3 



7 



(1 + 7)2 (1 + 67)2 

^ (174) 



max 

j = l,2.kd<K,l^k 



(1 + ^7)' 
1 







l|wo-P (1 + 7)^ 



(175) 
(176) 



as K,N ^ 00, § ^ p. 

Proof: The proof is similar to that of Lemma [3] the 
main difference lying in the fact that there is only one 
case to consider when looking at the interference terms (cf. 
( |135| l), since under coordinated beamforming, hj, ^ j is no 
longer independent of w^., j. Thus, defining j j = AI and 
Dj, J J = eAI, exactly as in the proof of Lemma [3] we can 
show that 



max 

j = l,2,k<K 



1 



trDfejj-A^ 



max 

j=l,2,k<K 



N 

— hfe,jjWfej - — trDfcjjAj 



> (177) 
0. (178) 



When studying the interference terms, the numerator in the 
RHS of ( |135| l, becomes, for the CBf case. 



N 



N 



'^k,j,i-Al,i,k^j,i^l i i^l,i,i-Ai,l,kJ,i^i^ 



(1 + Fhfcj,'A,,,,fc,j- ,h|^^. 

Considering the numerator in ( |179[ ), we can show that 



(179) 



max 

j,i=l,2,k,l<K,{k,j)^{lA) 



N 



'^kj,i-^l,i,kj,i'^l^i^i^l,i,i-^l,i,k,j,i'^k,j,i 



- ir7trDfej^,A,D,^,.iAi 



(180) 



Asymptotic analysis of the primal problem 

As in Appendix |Dj we start by assuming that the pkj 's in 
each cell are fixed at some common constant value pj, so that 
the total transmit power of cell j is Pj = f3pj. Once this 
regime is analyzed in the large system limit, we optimize the 
constants pi and p2 (to meet the DL SINR constraints) and 
show that the optimal constants are asymptotically optimal 
with respect to the primal optimization problem ( [T3| ). 

Under this regime in which = pj for all k in cell j, the 
following lemmas hold. 



We can also show that 



max 

j,i=l,2,k,l<K,{k,j)^{l,i) 



i^k,j,i-Ai,Lk,j,iii-k,j,i 



N 



trD 



kJA-l 



0. (181) 



The proof of the lemma is concluded by us ing the limits 
of the trace terms as derived in Appendix 
noting that for this specific case, all the 'Dkjj and D 



A, 
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matrices are equal, and are simply scaled identities, so that 
1 

rtrJJfcjjAj — 

eA 
N 
A 



N 
1 

N 
1 



N ' 



— trA, 



we therefore solve the following set of equations, 

1 _ li<?i^ 



7 



€7 



(182) 
(183) 



(1+7)" 



(1+67)2 



7 
e/3 



a.s. 1 



T7trA? 

7 



(1 + £7) 



2P2 + Pi 



P2- 



(1+7)" 



(l+e7)2 



A 1 + 



/3e 



(1+^7)' 



7 



1-/3 



(1+7)" 



7 
e/3 



1 A^ 
^trDfcj-,A,D,j,,A, - — trA^ 



Kit 



+7)" 



(1 + 67)^ 



(1+^7) 

(184) 



(185) 



(1 + £7) 



tPI + P2 



(i+7r 



(i+7r 



to obtain 



Pl=P2^P 
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(1+7) 



- ^7 
^l+£7 



This implies that 



TV 



-trD 



A2 



A,Di,,-,A, = — trA; 



(190) 



(191) 



(192) 



(193) 



£7 



(1+7)" 



(l+<;7)^ 



This choice of downlink transmit powers is guaranteed to meet 
(186) the SINR constraints as X, TV 00, % = 

We confirm the asymptotic optimality of this deterministic 
power allocation, together with the beamforming directions 
^ found from analysis of the dual problem, by verifying that the 
duality gap is tending to zero. Indeed, the primal objective 
Lemma 8. With Wkj = Akjhj^jj, such that A satisfies value converges to 2P ((j) = 1, since the power consumption 

and with pkj = Pj for k = 1, K, j = 1,2, it follows that, in both cells is the same), whereas the dual objective value 



with probability 1, 



E 



p. 



i',fe'.(fc',i')#(fcj) 



TV 



(1 + 7)2 (l + ejy 
(187) 



converges to 2/3X(t^. Recalling the definition of A, one can 
easily verify that these two quantities are the same, thereby 
completing the proof. 

J Equation ( |192| l shows that the primal problem has a solution 

/ 7 67 \ 
when effective bandwidth condition /? ( z — h 7 — I < 



Pj denotes the total transmit power of BS j, i.e. Pj — j3pj, 
J = 1,2. 

Proof: Since DL pkjS are all fixed to the same value, pj, 
we have 



E 



Pk'j' Ihfcjj'Wfc. 



j' ,k' ,{k' ,j')^{k,j) 



TV 



Pi 
N 



E 

k'=/^k 



llwfe'jIP 



TV 



^kj.j'^k'j I 



W 



(188) 



1+7 1 + £7 , 
1 is satisfied. Conversely, if this condition does not hold, there 

is asymptotically no feasible solution to the primal problem, 

in the limit as TV ^ 00. The latter observation follows from 

mono tonicity: The optimal powers are increasing functions 

of 7, but as the denominator of ( |192[ ) decreases to zero, the 

optimal total power from either BS, P, must diverge to infinity, 

and for higher values of 7 there can be no feasible solution. 

Appendix F 
Proof of Theorem[3] 

Throughout this section, let 



k'j\ 



Applying ( |176| ) and jllS) of Lemma |7] ( |188| ) becomes for 
fc = l,...,X, j = l,2. 



K 



A 



A 



K 



^EEh^.'h.'.'+I 



f=ik'=i 



Pj 



TV (1 + 7)2 
Noting that as K, TV -> cx 



+ e/3 



Pi 



oil). 



(189) 



E 

{k',j')Mk-j) 



K-1 

N 



(1 + 67)2 

/3 completes the proof. 



^kJ,Li 



E 

(k',j')^(k,j),{i,'i) 



^k',j'^k',f 



h^,,,,hk',,' 



(194) 



(195) 



(196) 



Referring to p6| ), and using results from Lemmas |7] and |8] 
to find the minimal pair of constants (pi,P2) for the two cells. 



where A will be specified later 

The following lemmas will be useful in proving the theorem. 
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Lemma 9. Let {^1,^2) satisfy < fJ-i < 2, i = 1,2, and 
define the function F{Xi,X2) by 

F(Ai,A2) = (Fi(Ai,A2),F2(Ai,A2)) 

/3 



^;(Ai,A2)=7 



1+7 

13 



(Aj + A -e) 



3 = 1,2. 
(197) 



r/ien 

i) F is an interference function 1127V . 

ii) // (3 — - — < 1 then there exists a unique solution 

'_ 1 + 7 

(Ai, A2) to the fixed point equation 

(Ai,A2) = F(Ai,A2). (198) 

Proof: Similar to the proof of Lemma [5] with A in ( |162[ ) 
replaced by 



1 + e 



1+7 



(199) 



Lemma 10. Given positive (/ii,/i2), <^nd (Ai,A2), define 
function 0(71,72) by 

G(7i,72) = (Gi(7i,72),G2(7i,72)) 

/3Aj e/3Aj 



Gj(7i,72) = 



Mi 



M.7 



1 



+ 7j 1 + Tj 



e/3A, 



1 



-7? 



1,2. 
(200) 



i) G is an interference function [27]. 

ii) If the inequalities 

71 > ^1(71,72), 72>G2(7i,72) (201) 

have a solution, then G has a unique fixed point. 

Proof: (i) is easily verified, and (ii) follows from lIZTl . 
Theorem 1. ■ 

Corollary 5. Let (^1,^2) satisfy < ^ii < 2, i — 1,2. 

7 ' ~ ~ 

Assume that (3- — — < 1, and let (Ai, A2) be the unique fixed 



1 +7 



point in ( |197|l , as identified in Lemma [9] TTzen the function G 
in Lemma [7Q| has a unique fixed point, namely (7,7). 

Asymptotic analysis of the dual problem 



Assume that (3 



7 



1+7 



< 1, let (/ii,/i2) be feasible for the 



dual (|29]l i.e. 111,112 > and /ii + /i2 = 2, and let (Ai, A2) 
be the unique solution to ( |198| l. Consider the suboptimal dual 
power vector that assigns power level Xj/N to all the users 
in cell j, j = 1,2, respectively. The dual UL SINR for any 
user k in cell j is given by Equation ( |34] i. The derivations 
leading up to ( |102| i in Appendix B-C show that with the given 



suboptimal dual power vector assignment ( [34] l will converge 
a.s. to (71,72) that is a fixed point of G as defined in 
Lemma 10 However, by Corollary |5] G has a unique fixed 
point, namely (7,7). It follows that the SINRs of all users in 
the system tend to the common value 7. 

We can follow the same reasoning as that in Appendix |E] 
replacing the optimization problem ( |165| l - ( |167| i by: 

max.^i^P2>o /3cr^ (Ai + A2) (202) 
s.t. ni+H2 = 2 (203) 
(Ai, A2) is the unique fixed point of ( |198| l. (204) 

The solutiorp^ to p02| i-( |204] i can be shown to be (/ii,/i2) = 
(1, 1), with corresponding (Ai, A2) = (A, A), where A is given 
in (|50|; The details of the proof are skipped due to space 
constraints. 

We thus conclude that the asymptotically optimal dual 
variables are /ii = /X2 = 1 and that all components of the 
optimal A corresponding to cell j users must converge to A, as 
given in ( [SO] ). Thus the optimal downlink beamforming vectors 
asymptotically point in the directions of the vectors given in 
( |53| ). Fixing the beamforming to lie in these directions, we 
now find the optimal power levels, pkj, to use for each user. 

Asymptotic analysis of the primal problem 

With the beamforming directions fixed, only the DL power 
levels pfej's still need to be determined. As in the proofs of 
Theorems [T] and |2] we first fix the DL powers to constants, and 
study the resulting limiting regime. Thus, we assume pkj 's to 
be fixed at a constant pj, j — 1, 2. Unlike the SCP and CBf 
cases, this does not imply that the total transmit power of BS 
j is equal to /3pj, but rather that the total transmit power of 
BS j is equal to 

Pj^^ llEjWfcjv 
N 



P. 



(205) 



"^^i^^i 11-^.' P 
Now turning to the DL power levels, we focus on the 
transmit strategy that allocates fixed power levels pj to all 
users in cell j. The following lemmas can be used to show 
that the left-hand side of ( (38] l converges in probability to cr^ 
with all users being allocated equal power as given by (|5TJ. 
This power allocation also results in zero duality gap. 



Lemma 11. With w 



Ak.jhj^ j, with X as given by ( |50| ." 



max 

j = l,2,k<K 



-(1 + e) 



1 - 



(i+7r 



a.s. 



(206) 



max 

j = l,2,k,k'<K,k'^k 



2e 
1+e 



1+7' 



^ P(l+e)2 (1+^)2 



-+0, i.p. (207) 



"Note that the function F in {i-91\ actually depends on the choice of 

(M1:M2). 
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max 

j = ia,k,k'<K,k'^k 



2e 



1 



Lemma 12. WiY/i w^j = Afe j-h^^, vv/f/; A as given by ( [SO] ), 
as K, N oo, ^ — > /3, f/ie following holds with probability 



1, for j = 1,2, 



^ 0, i> (208) 



2 ' 



as K,N ^ 00, § ^ l3. 

Proof: The proof is similar to that of Lem mas [3] and |7] 
Thus, defining Dfc.j and_ as in ( fT03] l and ( fT04| , ^1 ^'"fS^: 

= 1 and Ai A2 A, we can show, by applying Lemma Lemma [llj we can show that the ^||E^Wfcj|p and ^||wfe 
5.1 in 1 55 1 and Lemma [T] that 



(218) 



where Pj is the total transmit power under the considered DL 
transmit power regime (see \2Q5) ). 

Proof: Following similar steps to those in the proof of 



max 

j = l,2,k<K 



max 

j = l,2,k<K 



N 



\^kj 



P - ^trDz,,A2 



— h.k,j^ko 



1 

TV 
1 

TV 



trDfcjA 



Now consider the interference terms, 
{k' 7^ We have that 

2 



(209) 
(210) 
with 



converge a.s., and uniformly in their indices, to constants such 
that the following holds 

2 



max 

j,i=l.2,k<K 



0. 



(219) 



As a result. 



E 



2 



0(1). 



(220) 



\^k'3' 



A 



( 
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Similarly to ( |209| l and piO| l, we can show that 
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This completes the proof. ■ 
Thus, asymptotically, under this scheme, both BSs transmit 
with equal power. 

Lemma 13. With w^j = Akjhj^j, such that A satisfies jSQ) , 
and with pkj = pj for k = \, . . . ,K, j — 1,2, it follows that, 
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We can also show that 
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The proof of the lemma is concluded by us ing the limits 
of the trace terms as derived in Appendix B-C with Aj = A 
given by (|50|, 
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Proof: Since DL pfejS are all fixed to the same value, pj, 
we have 
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(215) Applying ^07) and ( p08] l of Lemma |7] (|222| for fc = 
1, . . . ,K, j = 1,2, will converge in probability to 
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(217) Noting that as if , iV 00, — ^ /3 completes the proof. 
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Referring to ( |38| ), and applying Lemmas [TT] and [13] to find 
the minimal pair of constants (pi,p2) for users in the two 
cells, we therefore solve the following set of equations, 
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to obtain 
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This implies that Pi and P2 (recall Lemma 12 1 both converge 
a.s. to P, where 

P:= f3p. (227) 

This choice of DL transmit powers meets the SINR constraints 
as K,N ^ 00, ^ — /3 with a probability tending to 1. 

We confirm the asymptotic optimality of this deterministic 
power allocation, together with the beamforming directions 
found from analysis of the dual problem, by verifying that the 
duality gap is tending to zero. Indeed, the primal objective 
value converges to 2P (0 = 1, since the power consumption 
in both cells is the same), whereas the dual objective value 
converges to 2/3 Xa"^. Recalling the definition of A, one can 
easily verify that these two quantities are the same, thereby 
completing the proof. 

Equation ( |226| l shows that the primal problem have a 
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solution when effective bandwidth condition (3 

is satisfied. Conversely, if this condition does not hold, there 
is asymptotically no feasible solution to the primal problem, 
in the limit as t 00. 

Appendix G 
Proof of Theorem[T] 



Using ( |43] l, r simplifies to: 
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which is positive at 7* = 0, and zero at 7* = 00. Defining 
f] = ^ + e and taking the derivative with respect to 7* we 

obtain - — equals 
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This has the same sign as 

/i(7*) = -log(l + 7*) + 



7*77 (1+7*) +7* 
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1+7^(1+7*)^ 

which is at 7* = 0, and —00 at 7* = 00. Differentiating the 
expression in (|230[), we obtain 
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This has the same sign as 

-(1 + ^(1 + 7*)')' 

+ (1 + 7*) (?? + 1 + ('/ + 1) 77 (1 + 7*)' - 27^7*2) 

= -7*\' - 7*^^ [1 + 377] - 7*^77 [1 + 377] - 7* [772 - 1] . 

If 77 > 1, then < V7* > 0, which implies hh*) < 
dj* 

V7* > and hence r(7*) is a decreasing function. If 7; < 1 
then h{j*) is positive for 7* small. This implies that r(7*) is 
increasing and then decreasing. 

Appendix H 
Proof of Theorem|2] 

The normalized achievable rate per cell is equal to 
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which is positive at 7* = and zero at 7* = 00. The sign of 
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which is at 7* = 0, and —00 at 7* = 00. The sign of 
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can be shown to be opposite to the sign of a polynomial in 
7* of degree 7, i.e. X]I=o '^^7**' ^^ch that 
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C3 = [2e (2e2 + 9e + 6) + l] + ae {e^ + 13e + 6) 

+ (l - - 16e) 
C2 = 3a^ [2e (e + 2) + 1] + a (Se^ + lOe + l) + e 

- (I7e^ + 8e^ + 4ae^) 
ci = (4e + 3) + 6ae + a - e (206^ + 3ae + Oe^ + 4e + 3) 
Co = a^ + 2ae-2ae2-2e3-e2_i^ (234) 

2 

where a = ^ . One can factor the constant term (independent 
of 7*) into 

(a + e + l)(a + e- 2e2 - 1). (235) 
We now consider two cases; 



A. a + e - 2e2 - 1 > 

Defining eo,i = ^^^^|^,eo,2 = ^+^|°~'^ , one can show 
that ( |235| ) will be positive if (a > 1, and < e < eo.2), or 
(| < a < li ™d £0,1 < e < £0,2)- One can further establish 
that if this is the case, all of the other coefficients of the 
polynomial ( |234| l will be positive and consequently /i(7*) 
will be a decreasing function, and hence /i(7*) < V7* > 0. 
In this case, r{j*) is a decreasing function and the optimum 
will occur at 7* — 0, which corresponds to /?* — 00. 

B. a + e - 2e2 - 1 < 

In this case, /i(7*) will be positive for 7* small. Thus, r(7*) 
will be increasing initially, but it decreases eventually (to zero), 
so it must attain its global maximum at a finite value of 7* > 0. 
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